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summary 


An  analytical  study  was  made  to  describe  the  structure  of  a strong 
normal  shock  wave  moving  into  argoi  and  the  nonequilibrium  flow  of  partially 
ionised  argon  in  the  laminar  boundary  layers  induced  behind  the  shock  wave 
on  the  shock- tube  walls  and  over  a flat  plate.  The  subsequent  interaction  of 
the  laminar  boundary  layer  with  a corner- expansion  flow  was  also  considered. 

In  order  to  determine  the  shock  structure,  the  ionization  and 

* f 

relaxation  processes  were  based  on  a two-step  model  of  the  collisional 
processes.  The  excitation  (or  ionization)  cross-sectiq/  constant  for  the 
argon  atom-atom  collisions  was'  determined  to  be  3*5  x (lC)  cr^/ev)  from  a 
comparison  of  theoretical  and  experimental  shock- structure  da  ;a. 

■ The  effects  of  a small  amount  of  hydrogen  inpurity  in  the  argon 
test  ,gas  on  shock- wave  structure  was  evaluated  and  discussed,  as  the 
hydrogen  inpurity  can  markedly  reduce  the  total  relaxation  length.  A study 
of  this-  effect  was  required  in  connection  with  stabilizing  the  experimental 
shock  waves , 

Aa  integral  method  was  used  in  the  analysis  to  study  both  the  shock 
induced  nonstabinnary  laminar  boundary  layer  on  the  shock- tube  walls  and  the 
quasi-steady  flat-plate  laminar  boundary  layer  for  an  ionized  argon  flow. 

The  frozen,  equilibrium  and  nonequilibrium  flow  solutions  were  obtained  and 
compared  with  some  existing  experimental  results.^ The  calculated  results 
based  on  the  integral  method  were  found  to  be  in  agreement  with  the  expert-  - 
mental  data  for  the  shock- induced  wall  boundary  layer,  but  they  did  not 
agree  with  the  experimental  data  for  the  flat-plate  quasi-steady  boundary 
layer . 

The  cold-wall  similarity  method  of  Sullivan  was  extended  to  apply 
to  the  interaction  of  a laminar  boundary  layer  with  a corner- expansion 
wave  for  a supersonic  frozen  flow.  The  validity  of  the  extended  method  is 
discussed. 
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Notation 


'fell 

n n polynomial  constant  for  the  velocity  profile  (Chapter  4) 

speed  of  sound  (Chapter  5) 

cross-sectional  area  of  a shock  tube 
fell 

n polynomial  constant  for  the  enthalpy  profile  (Chapter  -4) 

magnetic  field  intensity 

n n constant  used  in  Eq.  4 . 39"b 
fell 

n polynomial  constant  for  the  degree  of  ionization  profile 

(Chapter  4) 

ratio  P^/PgUg 
fell 

n constant  used. in  Eq..  4.39c 
defined  by  Eq..  5.2b 
specific  heat  at  constant  pressure 
defined  by  Eq..  3. 31 

skin  friction  coefficient  defined  by  Eg..  5.18 
ambipolar  diffusion  coefficient 
electron  charge 

1 

electric  field  intensity 

velocity  ratio  defined  by  Eq.-  4.22 

velocity  ratio  defined  in  Section  4.8 
till 

n integral  form  defined  by  Eq..  4 .23  (n  = 1,  2 or  3) 

total  enthalpy  ratio  defined  by  Eq.  4.22 

defined  by  Eq.  4.2 

defined  by  Eq.  5*3 

enthalpy  (specific) 

enthalpy  (total) 

ionization  energy 


•=25K^r3jra3aasi 


an  integral  defined  by  Eq.  4.44 

Boltzman  constant 

production  rate  used  in  Chapter  3 

equilibrium  constant  of  the  production  rate 

thermal  conductivity 

constant  defined  by  Eq.  4.23 

characteristic  length 

Lewis  number 

mass 

Mach  number 

particle  number  density 

pressure 

pressure  tensor 

Prandtl  number 

collision  cross-section 

u /u 
w e 

gas  constant 
Schmidt  number 

excitational  cross-section  constant 
time 

temperature 

particle  velocity  in  the  x-direction 
particle  velocity  in  the  y-direction 
drift  velocity 
mass  production  rate 

direction  parallel  to  the  shock- tube  wall 

direction  normal  to  the  shock- tube  wall 

ratio  of  the  degree  of  ionization  defined  by  Eq.4,22 


<2  degree  of  ionization 

£ defined  by  Eq.  4.50a 

7 isenfcropic  specific  heat  ratio  for  a partially  frozen  flow 

5 boundary  layer  thickness 

€ . internal  energy 

A defined  by  Eq,  4.27 

fi  viscosity 

v collisional  frequency 

p density 

H coordinate  defined  by  Eq.  4.20 

coordinate  defined  by  Eq.  4.62b 

£ coordinate  defined  by  Eq.  4.20 

!g  coordinate  defined  by  Eq.  4.62b 

potential  between  the  wall  and  the  plasma 


e 
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f 

H 

H+ 

i 

I 
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atom 
argon 
argon  ion 

electron  (Chapter  3) 

boundary  layer  edge  (Chapters  4 and  5) 

ionization  or  forward  reaction  (Chapter  3) 

hydrogen 

hydrogen  ion 

ion 

ionization 

recombination  or  backward  reaction 
sheath  region  (Chapter  4) 
shock  (Chapter  3) 
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1.  INTRODUCTION 


The.  structure  of  a strong 'shock  wave  in  an  ionizing  gas,  and  the 
associated  nonequilibrium  boundary  layer  flow  behind  such  a shock,  have 
received  considerable  attention  during  the  past  decade,  Research  into  shock 
structure  has  determined  how  elastic\and  inelastic  collisions  among  atoms, 
ions  and  electrons  can  influence  the  ionization  process*  Shock- structure 
measurements  present  an  opportunity  for  an  indirect  determination  of 
excitation  (or  ionization)  cross- sections  for  heavy  gas  particles. 

An  understanding  of  boundary- layer  flows* in  a partially  ionized 
gas  is  helpful  in  designing  a successful  reentry  into  'the  Earth’s  atmosphere 
at  supersonic  and  hypersonic  conditions.  If  the  term  boundary  layer  can  be 
applied  to  any  region  of  a continuum  flow  within  which  the (transport  of  mass, 
momentum  and  energy  by  diffusive  mechanises  is  .important,  then  the  boundary 
layer  of  an.  ionizing  gas  can  be  seen  to  be  generally  more  c duplex  than  that 
normally  encountered  in  a non-ionized  aerodynamic  type  of  flow.  The  presence 
of  ions  and  electrons  introduces  new  transport  mechanisms  in  the  boundary 
layer.  Additionally,  the  magnitude  of  the  various  transport,  properties  of 
a partially-ionized  gas  can  be  markedly  different  from  a perfect  gas  Even 
today,  after  many  years  of  research,  boundary  layer  flows  of  a partially-' 
ionized  gas  are  not  fully  understood,  experimentally  or  theoretically, 

Experiments  were  recently  conducted  at  UTIAS,  using  a hypervelocity 
shock  tube.  These  experiments  provided  uniqhe  and  reliable  data,  (inter- 
ferometric) on  both  the  structure  of  a strong  shock  wave  in  argon  and  on  the 
shock- induced,  nonequilibrium  boundary- layer- flows  in  ionized  argon,  in 
order  to  conpare  with  the  analyses. 

In  Chapter  2 the  'general  equations  of  motion  for  a partially 
ionized  gas  have  been  reviewed  briefly.  A simple  model  for  studying  initial 
ionization  and  relaxation  processes  behind  strong  shock  waves  in  an  argon 
flow  is  given  in  Chapter  3»  The  excitation  (or1  ionization)  cross-section 
constant  for  argon  atom-atom  collisions  can  be.  and  is  determined  from  a 
comparison  of  theoretical  and  experimental  results.  The  constant  determined 
herein  is  more  reliable  than  the  value  obtained  by  Kelly  (Ref.  12),  which 
has  been  widely  accepted  for  previous  theoretical  calculations.  The  im- 
purity effect  of  hydrogen  molecules  on  shock  structure  in  argon  is  also 
determined  and  discussed,  and  the  theoretical  and  experimental  results  are 
compared. 


The  laminar,  nonstationary,  shock  tube-wall  boundary  layer  behind 
a normal  shock  wave,  and  the  quasi- steady  fiat-plate  boundary  layer  in 
partially  ionized  argon,  are  both  considered,  in  Chapter  4.  The  method  of 
solution  was  based  on  the  Karman-Pohlhausen  integral  method.  The  integral 
and  similarity  solutions  are  compared  and  the  results  discussed.  Theoreti- 
cal calculations  for  nonequilibrium,  frozen  and  equilibrium  flows  are 
compared  with  UTIAS  experimental  data. 

In  Chapter  5 the  cold-wall  similarity  method  of  Sullivan  is 
adopted  and  extended,  in  order  to  treat  the  interaction  of  a laminar 
boundary  layer  with  a comer-expansion  wave  in  a partially  ionized  super- 
sonic argon  flow , The  flow  was  assumed  to  be  frozen  both  before  and 
after  the  corner  expansion.  Actually,  this  assumption  is  not  valid,  as 


deduced  from  experimental  results  (see  Chapter  5).  However,  the  assumption 
greatly  simplifies  the  calculations,  and  it  provides  a rough  estimation  of 
the  variations  of  pressure,  density,  temperature  and  Mach  number  as  a 
function  of  turning  angle  at  the  edge  of  the  boundary  layer.  The  validity 
of  this  oversinplified  model  is  discussed. 

In  the  last  chapter  discussions  are  given  on  the  significance  of 
the  theoretical  models. 


2.  GENERAL  EQUATIONS  OF  MOTION 

2.1  Introduction 

The  general  formulation  of  the  gasdynamic  conservation  equations 
for  individual  species  in  a nonequilibrium  partially  ionized  gas  mixture 
has  been  reviewed  by  Appleton  and  Bray  (Ref.  l) . The  conservation  equations 
for  the  electrons  and  the  overall  conservation  equations  were  developed  for 
a three- conponent  plasma  consisting  of  neutral  atoms,  singly  ionized  ions 
and  electrons.  In  this  formulation  the  main  assumption  is  that  each  of  the 
three  conponents  has  a Maxwellian  velocity  distribution.  Also,  the  electrons 
can  have  a tenperature  that  is  different  from  the  heavy  particles  and  they 
can  drift  relative  to  them. 

The  nonequilibrium  electron  tenperature  is  of  interest  in  gasdynamic 
problems  involving  ionized  gases,  namely,  ionization  and  relaxation  processes 
in  shock  waves  and  in  expansion  flows  in  rarefaction  waves,  nozzles  (Ref.  2) 
and  at  corners.  Nonequilibrium  phenomena  of  electron  (or  ion)  mass  production 
are  also  of  interest  in  shock  structure,  boundary  layer  and  expansion- wave 
(Ref.  3)  problems.  Recently,  Igra  (Ref.  4)  reviewed  briefly  the  relevant 
formulations  and  atomic  processes,  especially  the  three-body  recombination 
process.  The  latter  was  studied  in  some  detail. 

In  this  chapter  thb  basic  equations  are  presented  for  a nonequilibrium 
laminar  boundary- layer  flow  induced  behind  a normal  shock  wave  on  the  shock- 
tube  wall  or  over  a sharp  leading  edge  flat  plate.  The  basic  equations  for 
an  ionized  argon  plasma  flow  are  based  on  the  general  formulation  (Ref.  l) 
of  the  conservation  equations. 

2.2  General  Equations  of  Motion  for  an  Ionized  Gas 

An  ionized  monatomic  gas  or  plasma  is  considered  which  consists  of 
a mixture  of  atoms,  singly  ionized  ions  and  electrons.  For  each  Species  the 

Scrosccpic  balance  equations  canibe  expressed  (Ref.  l)  by  using  the  plasma 
croscopic  properties,  as  shown  below* 
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The  quantity  < 0 > is  the  average  of  the  property  0 , n is  the  number  density 
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velocity  of  a particle  of  species  s.  The  source  term  expresses  the  change 
in  < $s  > as  a result  of  both  external  influences  (i.e.,  electric,  magnetic 
and  gravitational  fields)  and  internal  influences  (i.e.,  chemical  reactions, 
heat  conduction,  diffusion  and  viscosity) . 

In  this  analysis,  for  a mixture  of  atoms,  ions  and  electrons,  it 
"Mill  he  assumed  that  each  species  has  a Maxwellian  velocity  distribution 
with  an  appropriate  temperature. 

The  equations  for  mass-production  rate,  momentum  and  energy  for 
the  electron  gas  are  given  below | 
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The  subscripts  e,  a,  and  i denote  electron,  atom  and  ion,  respect- 
ively j u and  Ug  are  the  velocities  of  the  heavy  particles-  (atoms  .and  ions.) 
and  electrons,  respectively;  we  is  the  rate  of  creation  (or  disappearance) 
of  electrons;  ma  and  mg  are  the  masses  of  the  heavy  particles  (mass  of  an 
atom  is  approximately  equal  to  that  of  an  ion)  and  electrons,  respectively; 
e is  the  electron  charge;  E and  B are  the  (electric  and  magnetic  fields, 
respectively;  e is  the  internal  energy;  vea  and  vej_  are  the  respective 
collisional  frequencies  between  electrons  and  atoms  and  between  electrons  and 
ions;  Qrad  is  radiation  source  term;  T and  Te  are  the  respective  tempera- 
tures of  heavy  particles  and  electrons;  pe  is  the  partial  (hydrostatic) 
pressure  of  electrons;  and  lx  is  the  net  energy  gained  by  the  electrons  per 
event  of  the  three-body  recombination  process. 

The  equations  of  continuity,  momentum  and  energy  for  the  whole 
plasma  are  obtained  by  summing  the  corresponding  equations  -for  all  the  plasma 
constituents.  The  following  definitions  and  relationships  are  used. 
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We  should  note  that  if  the  viscous  effects  can  he  neglected  (i.e.,  = 0) , 

then  Eqs.  2.8  and  2.9  reduce  to  the  equations  given  by  Appleton  and  Bray 
(Ref . l) . Additionally,  the  pertinent  set  of  equations  for  other  simpler 
problems  can  be  obtained  from  the  general  equations  - Eqs.  2.2  to  2.8. 

2.3  Equations  for  a One-Dimensional  Inviscid  Flow  of  an  Ionized  Gas 

For  studies  of  shock  structure  and  expansion  nozzle  flows,  the 
variation  of  dynamic  and  thermodynamic  quantities  in  the  direction  normal 
to  the  flow  direction  are  normally  small  compared  with  those  in  the  flow 
direction.  Hence  these  problems  can  be  treated  as  one-dimensional. 


Actually,  the  flow  for  the  shock-structure  problem  is  nonstationary. 
However,  a nonstationary  flow  can  be  readily  reduced  to  a steady  flow.  _Let 
(x,  y)  be  a coordinate  system  fixed  with  respect  to  the  wall,  and  let  (u,  v) 
be  velocities  parallel  to  (x,  y).  The  flow  is  unsteady  in  this  (x,  y)- 
coordinate  system.  Let  (x,  y)  represent  another  coordinate  system  which 


moves  with  the  same  speed  as  the  shock 
to  the  x-  and  y- coordinates  be  denoted 
coordinate  system  the  flow  is  steady, 
coordinate  systems  is  given  below. 

x = x - 

y = y 

u - u - 

V = V 


wave,  and  let  the  velocities  parallel 
by  u and  v respectively.  In  this 
The  transformation  relating  the  two 


u t 
e 

(2.11) 


The  velocity  of  the  wall,  uw,  equals  the  negative  value  of  the  shock  velocity, 
us.  Under  this  transformation,  nonstationary  flows  for  shock  structure  and 
boundary  layers  behind  a moving  shock  can  be  conveniently  treated  as  quasi- 
steady flows  (see  Fig.  l) . 


For  treating  the  shock- structure_prob^em_two  assumptions  are  made: 
(l)  no  electric  or  magnetic  fields  exist  (E  = B = 0),  (2)  ue  - u.  Note 
that  the  magnitude  of  the  electric  field  can  be  predicted  by  the  following 
expression. 


E1  a - 


n e 
e 


^Pe 

Sx1 


which  'is  negligible  for  the  present  study.  In  a similar  fashion  the  magnetic 
field  can  be  shown  to  be  even. less  important.  Due  to  the  very  small  electron 
mass,  the  effects  of  diffusion  (ue  u)  are  small  and  can  be  neglected.  The 
preceding  two  assumptions,  which  can  be  justified  for  many  types  of  nonequilib- 
rium shock- structure  problems,  have  been  widely  accepted  by  previous  researchers. 

For  the  present  shock- structure  study,  the  translational  transition 
region  of  the  shock  front  has  been  neglected,  since  its  thickness  is  negligibly 
small  compared  with  that  of  the  following  relaxation  "egion.  It  should  be 


noted  that  if  the  translational  transition  region  of  thesJ°ck  front  is 
considered,  gradients  in  the  flow  variables  are  large^and  the  effects^of 

viscosity  (p  ^)],  ^eat  conduction,  (Ka  3^)3  ( e 13x^' 


must  be  included. 

For  a steady  one-dimensional  and  inviscid  flow  of  a singly  ionized 
gas,  Eqs.  2.2  to  2.8  reduce  to'  the  following  expressions. 

Conservation  of  electron  mass: 


dn 


h (V1)  = at" 


dx 


(2.12) 


m 

<L  (1  n k T U)  = 3 n — (v  + v . ) k(T  - T ) 
^ ^2  e K Ae  5 e m-  ea  ei 


dn 

e 

dt 


dT 


dn 


I + <W  + k u[ne  dF  * Te  ST1 


(2.13) 


Conservation  of  total  mass: 


Total  momentum: 


dx 


(p  u A)  =0 


P u 


du  dp 

dx  dx 


Conservation  of  energy: 


dn 


p u fr  - u ll = - it  1 


(2.14) 

(2.15) 

(2.16) 


The  Boltzmann  constant  is  denoted  by  k,  A is  the  shock-tube  cross-sectional 
area  or  expansion  nozzle  cross-section,  I denotes  the  ionization  energy  of 
the  atom,  and 


P „ p£L  (T  +«Te), 


Hi  = (e  + p)/p, 


(2.17) 


€ = § [ne  k Te  + (na  + ne)  k T]. 


These  basic  equations  contain  five  dependent  variables:  %(x) , 

na(x) , T(x),  Te(x)  and  u(x).  Any  numerical  solution  will  depend  on  the  model 

adopted  to  describe  the  atomic  processes  for  and  Qrad* 

In  the  case  of  the  shock-structure  problem,  if  the  cross-sectional 
area  of  the  shock  tube  is  sufficiently  large,  then  the  effects  of  the  wall 
boundary  layer  on  reducing  the  total  area  will  be  negligibly  small.  > 
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A in  Eq.  2.14  can  be  taken  to  be  constant.  The  effect  of  the  wall  ‘boundary 
layer  can  be  included  simply  by  assuming  the  variation  of  the  cross-sectional 
area  outside  the  boundary  layer  is  a given  function  of  distance  (x) . Let 
Aq  be  the  initial  area  at  distance  x equal  to  zero,  then 


A(x)  = AqF(x)  (2.18) 

•where  F(x)  can  be  calculated  by  using  boundary- layer  theory. 

The  basic  difference  between  shockr structure  and  expansion  nozzle 
problems  is  in  the  atomic  collision  processes.  For  the  shock- structure 
problem  the  initial  ionization  processes  due  to  atom-atom  collisions  sire 
dominant.  By  contrast,  in  the  expansion-nozzle  problem,  the  three-body 
recombination  process  dominates  in  the  entire  relaxation  region.  The  former 
case  provides  an  opportunity  to  determine  the  excitation  .( or  ionization) 
cross-section  for  an  atom-atom  collision,  and  the  latter  case  allows  one  to 
determine  the  excitation  cross-section  for  electron-atom  collisions  or  the 
recombination  rate  constant. 

2.4  Equations  for  the  Laminar  Boundary  Layer  Flow  Induced  by  a Strong  Shock 
Wave 


For  a two-dimensional  problem,  Eqs.  2.2  to  2.8  contain  fourteen 
independent  variables:  ne,  na,  Te,  T,  u,  v,  ue,  ve,  Ex,  Ey,  Ez,  Bx,  By  and 

Bz.  Because  of  the  complexity  of  the  equations,  no  numerical  solutions 
appear  to  exist.  The  less  general  two-dimensional^  inviscid  plasma  flow,  with 
the  assumptions,  T = Te,  ue  = u,  ve  = v,  E = B = 0 and  Qraa  = 0 has  been 
solved  by  many  researchers.  For  example.  Glass  and  Takano  (Ref.  3)  deal 
with  the  problem  of  a corner- expansion  flow. 

The  occurrence  of  an  appreciable  degree  of  ionization  in  a gas 
flow  introduces  some  features  that  are  markedly  different  from  those  normally 
encountered  in  a perfect-gas  flow  or  in  a flow  with  chemical  dissociation, 
primarily  because  of  the  presence  of  ions  and  electrons.  The  reasons  are 
given  as  follows: 

1.  The  extremely  low  mass  of  the  electron  yields  a species  possessing  a 
thermal  conductivity  that  can  be  much  higher  than  that  of  any  other 
species  present  in  the  mixture. 

2.  The  collisional  energy- transfer  process  between  electrons  and  heavy- 
particles  is  relatively  slow  and  gives  rise  to  the  possible  situation 
in  which  the  electrons  may  have  a temperature  different  from  that  of 
the  heavy  species. 

The  detailed  analysis  of  an  appreciably  ionized  gas  is  necessarily  complex. 

In  general,  it  is  well  known  that  three  rather  distinct  regions 
exist  near  the  surface  of  a plasma  flow: 

1.  Away  from  the  wall  the  gas  is  quasi-neutral,  the  ion  diffusion  velocity 
is  small,  and  the  behaviour  of  the  ga.c  in  this  region  is  described  by 
continuum  equations. 
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2. 


Near  the  surface  hut  not  adjacent  to  it  the  gas  remains  quasi-neutral, 
hut  the  ion  diffusion  velocity  is  comparable  to  the  ion  sound  velocity 
and  the  usual  equations  are  not  valid  (the  region  will  he  referred  to 
as  the  "transition  region" . 

3.  Adjacent  to  the  surface  of  the  wall  a space  charge  sheath  exists  within 
the  gas  which  is  no  longer  quasi -neutral.  . • 


A major  difficulty  exists  in  connection  with  the  latter  two  regions,  because 
it  is  necessary  to  solve  the  Boltzmann  equation  for  each  species.  In  Fig.  2 
the  flow  regimes  near  a cold  surface  and  .the.  order  of  thickness  of  each  . 
regime  are  indicated.  ■ 


In  the  case  of  a partially  ionized  gas,  a sheath  formed  next  to  the 
wall  is  thin  compared  to  the  boundary  layer  thickness.  Hence,  the  ionized 
gas  in  the  boundary  layer  should  be  electrically  neutral  and  the  concept  of 
ambipolar  diffusion  can  be  adopted  (described  in  Section  4.3). 

The  present  section  is  specifically  concerned  vithi.the  basic 
equations  for  a nonequilibrium  boundary- layer  development  on  the  cold  wall 
surface.  The  basic  assumptions  are:  (l)  steady  flow,  (2)  laminar  flow 

(3)  no  continuum  radiation  losses,  (4)  all  species  have  the  same  mass 
motion  velocity,  (5)  free-stream  conditions  are  constant  along  the  flow 
direction,  (6)  no  electric  or  magnetic  fields  exist*  (7)  T = Te.  Note  that 
when  the  boundary  layer  is  inside  the  relaxation  zone,  assumption  (7)  may 
be  invalid.  Under  these  assumptions,  the  basic  equations  for  the  boundary 
layer  flow  are  given  below  (Ref.  47). 


Continuity  equation: 


Momentum  equation: 


is  /pu> + iy  (pV>  = 0 


du  , 8u  dp  , 5 ( „ 5u  'N 

pu3s  + pv^  = -d£  + ^VP^y 


Energy  equation: 


Sh  , Sh  5 
pu  3S  + pv  3S  = 


5 f u2 

" qc  " qa  + p <5y  \ T 


Conservation  of  species: 


5a  , 5a  5 r „ t , * 
PK  + pv  5?  * ^ t-  pi  V + wi 


(2.19) 


(2.20) 


(2.21) 


(2.22) 


In  these  expressions,  the  respective  symbols  p,  H,  qc,  q&,  V±,  Pi  and  w-^ 
denote  viscosity,  stagnation  enthalpy,  conductive  heat  flux,  diffusive 
energy  flux,  ion  diffusion  velocity  normal  to  the  wall,  ion  density,  and 
net  production  rate  of  ions. 


For  an  electrically  neutral,  singly  ionized,  monatomic  gas  consisting 
of  atoms,  ions  and  electrons,  all  at  the  same  temperature,  the  equilibrium- 
thermodynamic  relations  and  the  equation  of  state  for  the  mixture  are  given 
below. 
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p = p r T (l  + a) 


(2.23) 


c .,2 

H = | R T (l  + a)  + a I + — (2.24) 

where,  I denotes  the  ionization  energy.  Note,  the  enthalpy  due  to  excitational 
effects  is  neglected. 

The  conductive  heat  flux  (qc)  is  given  as  follows: 


q = - k 


(2.25) 


where  Cp  is  the  "frozen"  specific  heat  at  constant  pressure  and  defined  below. 


cp  = | R (1  + os) 


(2.26) 


R is  the  gas  constant  for  the  atom  and  k is  the  thermal  conductivity  for  the 
mixture  if  no  chemical  reactions  took  place. 

The  diffusive  energy  flux  for  ambipolar  diffusion  (q^)  is  given  by 
the  following  expression:  1 


9*  = - P [I  + i B T] 


(2.27) 


where  is  the  ambipolar  diffusion  coefficient. 

The  diffusive  mass  flux  of  ions  (Pi  Vi)  for  ambipolar  diffusion  by 
ion-electron  pairs  is  governed  by  Rick's  law,  as  given  below. 


pi  Vi  = " p D 


da 
am -By 


(2.28) 


By  introducing  the  Prandtl  and  Lewis  numbers  defined  by 


r ‘ K 


(2.29) 


p D c 
am  t 


the  ion  (or  electron)  concentration  and  energy  conservation  equations  take 
the  following  form: 


da  J da  d [ ’ Le  da  1 • 

pu3x  + pvBy  = By[^P'Byj  +wi 


(2.30) 
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(2.31) 


pu  1e  + pv  sf 


5 [n  5h  (.  i \ b ( t 


- l)  (I  + | R T) 


ba 

By 


In  general  the  Prandtl  and  Lewis  numbers  are  functions  of  the  gas  temperature 
(T)  and  degree  of  ionization  (a),  as  shown  in  Section  4.2.  This  temperature 
and  degree  of  ionization  dependence  couples  Eqs.  2.19  and  2.20  to  Eqs.  2.30 
and  2.31,  making  it  difficult  to  obtain  a solution. 


3 . SHOCK  STRUCTURE  AND  INITIAL  IONIZATION  PROCESSES 


3*3.  Xax  production 

Experimental  and  theoretical  investigations  of  ionization  rates  and 
relaxation  processes  behind  strong  shock  waves  moving  in  monatomic  gases  can 
be  found  in  Refs.  5 to  19.  It  is  now  well  established  that  when  a strong 
shock  wave  raises  the  temperature  of  atoms  from  a low  temperature  (about  300*K) 
to  a high  temperature  (over  10,000°K) , the  initial  ionization  is  due  to  atom- 
atom  collisions.  The  rate  of  ionization  is  controlled  by  the  rate  of  excitation 
from  the  ground  state  to  the  first  excited  level.  Ultimately  when  the  number 
of  electrons  becomes  sufficiently  large,  electron-atom  and  ion- electron- electron 
collisions  dominate  the  excitation,  ionization  and  recombination  processes  (see 
Eq.  3-1)  j and  the  rate  of  ionization  is  controlled  by  the  rate  of  energy  trans- 
fer between  the  heavy  particles  and  the  electrons.  The  radiative  processes 
will  be  important  (Ref.  5)  for  a low  electron  number  density  (e.g.,  at  T rr 
3200°K  and  ne  < 1012  cm-3) . Treatments  of  relaxation  phenomena  in  radiating 
argon-plasma  flows  can  be  found  in  Refs.  6,  7>  8 and  9.  For  specific  conditions 
of  a shock  Mach  number  of  15  and  an  initial  pressure  of  1 torr,  Kamimobo  et  al 
(Ref.  9)  have  shown  that  the  effect  of  radiation  on  the  relaxation  profiles  is 
negligible.  Oettinger  and  Bershader  (Ref.  6)  have  shown  that  the  effect  of 
radiation  is  negligibly  small  only  until  the  end  of  the  relaxation  zone. 
Thereafter  radiative  emission  becomes  appreciable. 

In  general,  the  excitation  cross-section  for  atom-atom  collisions 
is  well  known.  The  shock  tube  presents  an  opportunity*  for  indirect  measure- 
ments to  be  made  of  the  excitation  (or  ionization)  cross-section  of  heavy  gas 
particles.  Recently,  Harwell  and  Jahn  (Ref.  10)  have  employed  a transverse 
microwave  probe  bo  determine  the  cross-section  constant  of  proportionality, 

Sg,  as  equal  to  7 x 10~19  cm2/eV,  for  argon  inelastic  atom-atom  colli- 
sions. Morgan  and  Morrison  (Ref.  11)  have  made  a theoretical  reassessment 
of  the  ionization  mechanism  and  referred  bo  earlier  experimental  measurements. 
They  showed  that  a best  fit  curve  to  the  ionization  relaxation  time  measure- 
ments of  Petschek  and  Byron  (Ref.  5)  was  obtained  by  reducing  Harwell  and 
Jahn's  value  (Ref.  10)  for  S^r_^r  by  a factor  of  ten.  Kelly  (Ref.  12)  has’ 
reduced  the  impurity  level  for  his  experiments  and  obtained  a value  for  . 

£>Ar-Ar  of  1.2  x 10  ± 15$  cm2/eV,  which  has  been  widely  accepted  by  many 

researchers  (Refs.  2,  7>  10  and  13)  for  their  theoretical  calculations. 

However,  we  must  point  out  that  there  are  two  important  effects  which  must 
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be  considered  in  regard  to  Kelly's  experiments:  (l)  the  range  of  the  shock 

Mach  number  was  between  7 and  10,  which  is  too  low,  and  (2)  the  relatively 
small  cross-sectional  area  of  the  shock  tube  (5cm  x 5cm)  makes  the  wall 
boundary  layer  effect  important  (i.e.,  dA/dx  j=-  0) . On  the  other  hand, 
McLaren  and  Hobson  (Ref.  i4)  have  used  double  electrostatic  probes  to 
measure  t^  ionization  rates  and  they  obtained  a value  for  of 

2.5  x 10"  cm^/eV,  which  is  lower  than  the  value  suggested  by  Morgan  and 
Morrison  (Ref.  11).  Although  many  experiments  have  been  made  with  argon 
to  determine  the  initial  ionization  processes,  the  excitation  cross-section 
for  argon  atom-atom  collisions  has  not  been  determined  with  great  accuracy. 

An  experimental  investigation  (Ref.  15)  of  the  initial  ionization 
processes  in  a strong  shock  wave  in  argon  and  in  a mixture  of  argon  and 
hydrogen  (0.4$  by  partial  pressure)  was  mad,e  recently  in  the  UTIAS  4-in 
x 7-in  hypervelocity  shock  tube.  This  investigation  provided  the  most 
reliable  data  for  determining  Sjtr-Ar*  Considerable  care  was  taken  to 
reduce  the  impurity  levels  to  approximately  10-4  torr.  Therefore  at  an 
initial  pressure  of  about  5 torr  for  the  present  experiments  the  impurity 
level  was  about  20  parts  per  million  or  less.  The  large  cross-sectional 
area  made  the  effect  of  the  wall  boundary  layer  negligible  during  the 
experiments.  The  purpose  of  the  present  chapter  is  to  propose  a simple  j 
theoretical  model  from  which  the  excitation  cross-section  constant  of 
proportionality  for  argon  atom-atom  collisions  can  be  determined  from  the 
experimental  results  (Ref.  15) . The  effect  of  the  small  hydrogen  impurity 
in  the  argon  test  gas  on  the  ionization  rate  is  calculated,  compared  with 
experimental  work  (Ref.  15)  and  discussed. 


In  the  present  study  the  transition  through  the  translational 
shock  front  is  neglected  since  its  thickness  is  negligibly  small  compared 
with  that  of  the  relaxation  region.  In  the  translation  transition  regior 

d ( „ du 


in  the  flow  variables  are  large  and  the 


and  heat  conduction 


.effects  of  viscosity 


. Ie  (*a  Ie)  b ( Ke  |l  / 


must  be  considered. 


However,  as  only  the  relaxation  region  is  considered,  where  the  gradients 
are  small  (see  Figs.  6 and  7) , the  viscous  and  heat  conduction  dissipative 
processes  can  be  neglected.  Then  the  only  important  processes  in  the  flow 
are  the  collisional  processes-. 


3 .2  Collisional  Ionization  Rates 


Petschek  and  Byron  (Ref.  5)  and  Harwell  and  Jahn  (Ref.  10)  have 
shown  that  excitation  from  the  ground  state  to  the  first-excited  state  is 
rate  controlling  for  the  overall  ionization  process.  This  two-step  process 
is  based  on  the  fact  that  the  cross-section  for  excitation  from  the  ground 
state  is  greater  than  that  for  ionization  from  the  ground  state.  This 
familiar  two-step  collisional  ionization  model  is  adopted  herein.  Of  course, 
the  multi-step  collisional  radiative  ionization  models  of  Hollenbach  and 
Salpeter  (Ref.  1 6),  Kamimoto  et  al  (Ref.  9)  and  Bates  et  al  (Ref.  17)  are 
the  most  accurate.  However,  many  unknowns  are  contained  in  those  models 
and  some  simplifying  assumptions  are  required  before  actual  calculations 
can  be  made.  For  example,  there  are  three  unknowns  in  the  ladder- climbing 
model  of  Hollenbach  and  Salpeter,  where  transitions  to  neighbouring  levels 
are  only  allowed.  Kamimoto  et  al  have  shown  that  the  results  for  argon 
atom  and  electron  number  densities  as  calculated  by  a two-step  miodel  are 
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nearly  the  same  as  those  calculated  using  a multi-step  model,  except  that 
the  relaxation  time  based  on  the  multi-step  model  was  somewhat  longer. 

The  collisional- ionization  kinetics  of  singly  ionized  argon  are 
postulated  to  take  place  according  to  the  following  equations: 


Ar 

+ X 

- 

A*  + X 

(3.1a) 

A* 

+ X 

C“ 

A+  + e + X 

(3.1b) 

r 

r 

Ar 

+ X 

kf> 

kT 

A+  + e + X 
r 

(3-lc) 

The  two-step  model  is  represented  by  Eqs.  3 .la  and-  3* lb  for  the  overall 
reaction  path  given  by  Eq.  3.1c;  X denotes  either  the  argon  atom  (Ar)  or  an 
electron  (e). 

For  this  collision  process  one  can  express  the  electron-production 
rates  (iL)a  and  (ne)e  in  terms  of  recombination  rates  and  equilibrium  constants 
(Ref.  13),  as  illustrated  below: 

«.•>.  ■ nl  - kra  \ ne 

= kra  ^a)  na  ^Keq 

= kfe  na  ne  ' kre  ne 
* kre  <Te>  "e  <*.> 


n - n?]  (3.2) 

a e 


n - n®]  (3.3) 

a eJ 


The  electron-production  rate  due  to  atom-atom  collisions .is  denoted  by  (ne)a 
and  that- due  to  atom-electron  collisions  is  denoted  by  (ne).e,  and  the 
subscripts  a and  e denote  atoms  and  electrons,  respectively.  The  two  equi- 
librium constants  are  defined  by  the  following  expressions: 


K (T  ) 
eq  ' a' 


nt  „ (Tj 


n 


2 

e,eq  '~a,' 

(T  ) 
a,eq  x a.' 


(3.4a) 


(Te) 


ns  (T  ) 
e,eq  x e/ 

ne,eq 


(3.4b) 


The  calculations  of  kya  and  kre  require  a knowledge  of  the  depend- 
ence of  the  inelastic-collision  cross-section  for  the  first-excited  level 

(o*  ) and  on  the  kinetic  energy  (ej . A reasonably  good  approximation  is 
- ax'  , , x 

given  below: 


(3.5) 


a* 

ax 


S* 

ax 


(e  - €*)  with  6 > e* 
' x ay  x a 


where  Sax  fc^e  constant  of  proportionality  of  the  first-excitation  collision 
cross-section  between  atom  a,  and  particle  x,  and  ea  is 'the  excitation  energy 
of  the  first  level. 


By  applying  Eq.  3*5,  kra  and  kre 
and  13)  in  terms  of  for  argon  atoms  as 


can  be  expressed  (as  in  Refs, 
shown  below: 
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kra(Ta)  = 4.835  x 10 


.18 


S* 

°Ar-Ar 


T* 

T" 


+ 2 


(3.6) 


k (T  ) = 1.843  x 10"15  s* 
re'  e'  ^ Ar-e 


exp 


_„6 

cm 

sec 


(3.7) 


T*  and  Tp  are  the  first- excitation  and  ionization  temperatures  respectively 
for  an  argon  atom,  and  S*  is  in  units  of  cmP/eV.  (Rote  that  kra  in  Eq,  3.6 
has  now  been  divided  by  a factor  of  2 to  avoid  counting  like- like  collisions 
twice.) 


A schematic  diagram  of  the  collision  processes  from  the  Rankine- 
Hugoniot  translational  shock  front  through  the  relaxation  region  appears  in 
Fig.  3.  The  boundary  layer  formation  in  this  region  is  also  indicated. 
Radiation  losses  give  rise  to  a somewhat  nonuniform  equilibrium  flow.  It 
is  important  to  take  this  loss  into  account  during  the  calculation  of  the 
physical  properties  of  the  flow. 

3*3  Governing  Equations 


The  theoretical  approach  to  solving  the  shock  structure  is 
similar  to  that  described  in  Ref.  13,  except  the  correct  equations,  given 
in  Section  2.3  are  used.  Rote  that  the  conservation  equations  of  electron 
energy  described  in  Ref.  13  are  only  an  approximation  for  Eq.  2.13.  The 
governing  equations  are  obtained  from  Eqs.  2.12  to  2.l6,  as  given  below: 


dn 

-r~  (n  u)  = -r— 
dx  ' e ' dt 


= n 


(3.8) 


v + v . 
ea  ei 

u 


2T  . 

(T  - T ) - — ^ 

V e;  3u  dx 


un 


(t 


+ i V 


(3.9) 


pu  = PxUx 


(3*10) 
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(3.11) 


p + ,pu2  = pi  + PlUi 


| R(T  + o®  ),  + C®Tt  + ~ a2  = | RTi  + ^ 


(3-12) 


The  subscript.  1 refers  to  the  state  of  the  gas  in  front  of  the  translational 
shock  front,  vea  and  ve^  denote  the  elastic  collision  frequencies  for  the 
pairs  electron- atom  and  electron-ion  respectively,  which  can  be  expressed 
in  terms  of  the  elastic-ccillision  cross-sections  crel  and  c j^I,  as  shown  below: 


ea 


ex- 


v = U-OpP 

ea  m 

a 


/ 8kT  \k 
( - j a 

\™  J 


.el 

ea 


(3.13a) 


V . = — 

ei  m V 77m  J ex 


a s e 

The  degree  of  ionization  (a)  is  defined  by  Eq.  2.5. 


(3.13b) 


From  the  fact  that  ne  = ne(d/a),  Eqs.  3.2  and  3*3  yield  the  following 
expressions,* 

r of  (T) 

(1  - a)  - a2  (3.14a) 


da(T)  = (1  - a) 


m 


a J 


kra(T) 


eq' 


“e<Te>  - “ ( ? J kre<Te>  [ ll^'TfT  (l  ’ “>  ’ 
a y «-  eq  e'  1 

The  conservation  of  electron  mass,  Eq,  3.8,  becomes 


cr 


(3.14b) 


~ = (a  + d )/u 

dx  ' a e" 


(3.15) 


where  aa  and  (Xe  denote  production  rates  due  to  atom-atom  and  atom-electron 
collisions,  respectively. 

The  equilibrium  value  of  the  degree  of  ionization,  ae^(T) , can  be 
(.jtained  from  the  following  equation  (Ref.  3): 


rrsb  - r [ h T exp(‘  VT> 


cr  p, 

"eq 


(3.16) 


The  characteristic  density  for  ionization  pj  is  defined  by  Eq.  2.2^  of  Ref.  3 
(see  Table  1 of  this  report). 

*Note  that  Eqs . 15  and**l6  of  Ref.  13  are  incorrect. 
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The  value  of  du/dx  which  appears  in  Eq.  3*9  is  obtained  easily  from 
Eqs.  3*10  to  3.12,  and  it  is  given  below, 


du 

dx 


6uTi 

5 # Tl 
s 


>2  rp 

1_  A + 96a  _I  2 da 

M2  ' 5M2  T J 

s s 


(3.17) 


3.4  Model  Calculations  and  Comparison  with  Experimental  Results  for  Pure 
Argon 

For  the  present  calculations  the  same  collisional  cross-section 
employed  by  Petschek  and  Byron  (Pef.  5)  with  s|j._e  equal  to  7.0  x 10  18  cm2/eV 
was  used.  The  values  of  cr|*  have  been  taken  from  Refs.  11  and  13.  The  initial 
conditions  employed  for  the  calculations  for  a shock  wave  in  pure  argon  are 
summarized  in  Table  2,  which,  coincide  with  those  used  in  the  ’experimental  work 
in  the  UTIAS  4-in  x 7-in  hypervelocity  shock  tube  (Refs.  15  and  18) . 

The  first-order  differential  equations  (Eqs.  3*9  and  3.15)  were 
solved  by  using  a standard  Runge-Kutta  method  (see  Appendix  A) . The  initial 
condition  for  a was  taken  as  a(0)  a 0.  The  initial  value  of  the  electron 
temperature  is  somewhat  ambiguous.  Three  values  of  the  initial  electron 
temperature  are  given  as  follows: 

1.  Te(0)  = Ti(~  300°K), 

2.  T (0)  = T(0),  and 

3.  local  steady- state  value. 


However,  it  should  be  noted  that  these  different  initial  values  would  not 
affect  the  values  of  the  other  physical  quantities,  as  the  initial  number 
density  of  the  electrons  is  very  small,  or  a is  approximately  zero. 

For  a strong  shock  wave  moving  in  pure  argon  (Case  I of  Table  2), 
predicted  and  experimental  variations  of  the  degree  of  ionization  (a)  with 
distance  (x)  through  the  relaxation  region  are  shown  in  Fig.  4.  The  degree 
of  ionization  initially  increases  rather  slowly  from  zero  at  the  shock 
front  (x  = 0)  and  eventually  rises  rapidly  to  its  equilibrium  value  (ae  = 
0.158) . The  dashed  line  corresponds  to  a solution  using  the  value  of 
sAr-Ar  elUal  t 0 °«6  x 10" 19  cm2/eV.  It  can  be  seen  that  this  predicted 
variation  for  a has  the  same  features  as  displayed  by  the  experimental  data 
(Ref.  15),  but  the  predicted  relaxation  length  (or  time^  is  significantly 
shorter.  It  would  be  even  worse  for  larger  values  of  S Ar-Ar*  addition, 

a number  of  calculations  with  different  S fr-Ar  values  showed  that  the 
relaxation  length  increased  with  decreasing  values  of  sj^.  /w..  The  continuous 
line  in  Fig.  4,  corresponding  to  a solution  with  a lower  S^,-Ar  value  of 
3.5  x 10  2°  cms/eV,  which  is  in  good  agreement  with  the  experimental  data. 
Hence,  from  such  a comparison  the  excitations!  cross-section  constant  for 
argon  atom-atom  collisions  (S^.^j.)  has  been  determined  to  be  3.5  x 10  2 
cnr/eV.  This  newly  determined  value  is  used  for  the  remainder  of  the 
calculations . 
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The  predicted  and  measured  variation  of  density  through  the  relaxation 
region  is  shown  in  Fig.  5*  The  predicted  and  measured  density  profiles  are  in 
good  agreement,  adding  further  verification  to  the  acceptance  of  the  new  value 
for 

The  atom- temperature  (T)  and  electron- temperature  (Te)  variations 
through  the  relaxation  zone  are  displayed  in  Fig.  6.  Although  three  different 
initial  electron  temperatures  were  selected,  the  electron- temperature  profile 
is  affected  by  this  choice  only  in  the  small  initial  part  of  the  relaxation 
region,  as  shown  in  the  figure. 

The  flow  velocity  and  pressure  variations  in  the  relaxation  region 
are  given  in  Fig.  7>  and  in  Fig.  8 the  variations  of  electron  number  density 
(ne)  and  ionization  production  rates  are  shown.  It  can  be  seen  that  the 
ionization  production  rate  due  to  atom- atom  collisions  (aa)  is  very  large  in 
the  initial  stage  of  the  ionization  process  as  compared  with  that  due  to 
atom-electron  collisions  (de) , and  the  collision  process  is  controlled  by 
atom-atom  collisions.  The  ionization  production  rate  de  increases  very  # 
rapidly  as  the  number  of  electrons  increase.  After  de  exceeds  aa,  then  da 
quickly  decreases  and  the  collision  process  is  controlled  by  atom- electron 
collisions.  Eventually  ae  reaches  a maximum  and  then  falls  off  rapidly,  and 
the  collision  process  is  controlled  by  three-body  recombination  (electron- 
electron-ion)  and  radiation  processes. 

Additional  predicted  and  measured  results  (Ref.  15)  for  the  density 
and  degree  of  ionization  variations  through  the  relaxation  zone  are  given  in 
Fig.  9.  These  results  for  Case  II,  Table  2,  are  not  too  different  from  those 
of  Case  I,  and  the  predicted  and  measured  data  are  in  faij.  agreement. 

For  Case  III,  Table  2,  for  which  the  shock  strength  and  degree  of 
ionization  are  both  less,  the  predicted  and  measured  degree  of  ionization 
profiles  for  the  relaxation  region  are  compared  in  Fig.  10.  The  agreement  is 
not  as  good  as  for  Cases  I and  II.  This  disagreement,  however,  is  most 
likely  due  to  inaccurate  experimental  measurements.  Brimelow  (Ref.  15)  has 
mentioned  that  it  was  difficult  to  obtain  accurate  data  from  interferograms 
that  exhibit  small  fringe  shifts,  as  is  the  case  when  the  degree  of  ionization 
is  small.  However,  it  should  be  noted  that  the  analysis  predicts  the  correct 
relaxation  length. 

t should  be  pointed  out  that  the  analysis  overpredicts  the  relaxa- 
tion lengths  as  measured  by  Oettinger  and  Bershader  (Ref,  6).  The  initial 
conditions  for  both  Brimelow' s work  (Ref.  15)  and  Oettinger  and  Bershader 
(Ref.  6)  were  the  same.  However,  the  measured  relaxation  lengths  of  Oettinger 
and  Bershader  were  shorter.  The  reason  for  the  different  lengths  is  not  clear. 
The  shock  tube  used  by  Oettinger  and  Bershader  had  a smaller  cross-section  and 
its  impurity  level  is  not  known. 

In  some  of  the  interferometric  results  of  Brimelow  (Ref.  15),  it 
was  found  that  the  degree  of  ionization  was  larger  near  the  shock- tube  wall 
than  at  the  centre  of  the  frees tream,  and  relaxation  lengths  were  shorter 
near  the  wall.  At  the  present  time  it  is  difficult  to  explain  this  phenomenon. 
It  may  be  possible  that  a two-dimensional  model  incorporating  the  viscous 
effect  and  impact  phenomena  between  ions  (or  electrons)  and  a solid  surface 
would  be  required  to  explain  the  above-mentioned  phenomenon,  or  perhaps 
impurity  gradients  (e.g.,  water  molecules)  closer  to  the  wall  may  account  for 
it. 
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3.5  Effects  of  Hydrogen  Impurity  on  the  Ionizafcion  Processes 

The  effects  of  inpurities  in  the  test  gas  on  the  structure  of  shock 
waves  was  studied  initially  by  Morgan  and  Morrison  (Ref.  11) . In  this  section, 
the  effects  of  a small  addition  of  hydrogen  (impurity)  to  the  argon,  test  gas 
on  the  shock  wave  structure  is  considered  in  some  derail. 

The  initial  ionization  of  hydrogen  in  a strong  shock  wave  has  been 
studied  by  Belozerov  and  Measures  (Ref.  19) . From  a comparison  of  theoretical 
and  experimental  results  they  have  determined  to  be  3*57  x 10-18  cm2/eV, 
which  is  about  l/l5  that  of  the  corresponding  cross- section  for  electron-atom 
excitation  collisions  where  S^_e  equals  5*1  x 10” 17  cm2/eV.  They  showed  that 
the  electron  temperature  is  nearly  identical  to  the  atom  terrperature  every- 
where in  the  relaxation  region  except  for  a very  small  region  near  the  shock 
front  (x  ~ 0) . This  result  is  quite  different  from  a similar  solution  for  an 
argon  test  gas  (see  Section  3.4).  The  reasons  are  as  follows:  first,  the 
mass  of  the  hydrogen  atom  is  l/40  that  of  the  argon  atom,  and  second,  there 
is  no  Ramsauer  effect  (Ref.  19)  for  the  electron-hydrogen  atom  elastic  colli- 
sions. The  latter  effect  reduces  the  cross-section  for  elastic  electron- argon 
atom  collisions. 

Chang  (Ref.  20)  has  shown  that  for  certain  shock  velocities  it  is 
possible  to  consider  the  different  relaxation  processes  (e.g.,  translation, 
vibration,  dissociation,  and  ionization)  independently  for  hydrogen  molecules. 
For  the  present  work  we  can  assume  that  ionization  of  the  hydrogen  molecules 
behind  a strong  shock  starts  effectively  after  dissociation  is  coirplete. 

Also,  the  relaxation  length  for  the  dissociation  process  is  very  small 
compared  with  ionization.  With  these  assumptions  there  appears  to  be  no 
difference  between  hydrogen  and  monatomic  argon  in  the  theoretical  description 
of  the  ionization  process . The  only  difference  in  the  ionization  processes 
between  pure  hydrogen  and  argon  is  that  the  initial  conditions  for  hydrogen 
will  correspond  to  those  for  a fully  dissociated  gas  behind  the  shock  front. 
These  required  conditions  have  been  given  by  Belozerov  and  Measures  (Ref.  19) . 
For  the  present  case  (Case  IV  of  Table  2),  since  we  will  consider  only  a 
small  amount  of  hydrogen  impurity  (0.4%  by  pressure)  in  the  argon  test  gas, 
the  assumption  is  made  that  the  initial  conditions  for  hydrogen  and  argon  are 
identical. 


The  following  reactions  are  considered  for  the  collisional  ionization 


processes . 
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Processes  (a)  and  (b)  for  pure  argon  have  been  discussed  previously  in  Section 

3.4. 

The  production  rates  due  to  atom-atom  collisions  (a,  c,  d,  and  e) 
and  atom-electron  collisions  (b  and  f)  can  be  expressed  as  follows: 
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where  it  should  be  noted  that  (n  )Ar-H  ^ (A  )H-Ar, 

' e a N e a 

A B A-  6 

The  forward  rate  coefficients  kfa  between  atoms  A and  B and  kfe 

between  atom  A and  an  electron  can  be  written  in  terms  of  the  excitational 

cross-section  constants  S and  sj£_e  as  shown  below  (Ref.,  11): 


4'z  <T>  = SI-B 


3Z(V*B\ 

“A 


m.  <./  _ 


1/2  3/2  / T*  \ / -TJ 

(kT)  f 2T  + 1 ) sxp  ( T 


(3.20) 
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(3.21) 


These  rates  must  be  divided  by  two  for  like-like  collisions  to  avoid  a double 
count . 

The  forward  and  recombination  rates  are  related  to  the  equilibrium 
constants  as  indicated  below; 


kA_B  = kA-B  (T )/KA_B  (T) 
ra  fa  ' " eq  v ' 


(3.22) 


kr'e6  = 4?  (Ta><e  <Te> 


(3-23) 


'■’wr^owi 


The  equilibrium  constants  are  defined  below: 
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follows: 
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e»  e 

The  definition  of  the  degree  of  ionization  for  the  mixture  is  as 


ne,Ar  + ne,H 
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If  we  define  a ratio  C*  as  shown  below, 
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then  a can  be  expressed  in  terms  of  C*  as  follows: 
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where  and  are  defined  below, 
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For  the  present  calculations  the  excitation  collision  cross-section  fbr  atomic 
hydrogen  by  electron  impact  is  given  by, 


and 


°H-e  = 5,1  x 10"17  (10,2  “ E^CI"2 


(3.35) 
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as  obtained  by  Belozerov  and  Measures  (Ref.  19)  from  a comparison  of  theoretical 
and  experimental  results.  S^r_pj  and  are  still  unknowns  for  the  present 

work.  However,  the  method  used  by  Kelly  (Ref,  12)  in  his  treatment  of  the 
argon-xenon  case  by  assuming  = S|r_Ar  and  = S^_H  is  adopted  herein. 

This  assumption  would  be  invalid  for  a high  impurity  level  of  H due  to  the  small 
mass  ratio  between  H and  Ar.  However,  it  can  be  accepted  for  the  present  work 
as  the  impurity  level  in  the  test  gas  (argon)  is  small. 

Calculations  were  made  for  Case  IV  of  Table  2 for  the  relaxation 
zone  behind  a strong  shock  moving  in  argon  with  a small  amount  of  hydrogen 
(0.4%  by  pressure).  The  resulting  density  and  degree  of  ionization  variations 
through  the  relaxation  zone  are  given  in  Fig.  11,  along  with  the  experimental 
data  (Ref.  15).  The  agreement  between  the  predicted  and  measured  results  is 
good,  confirming  the  choice  of  the  theoretical  model.  From  a comparison  of  the 
results  of  Figs.  9 and  11,  it  is  readily  seen  that  the  total  relaxation  length 
of  the  pure  argon  shock  transition  is  substantially  reduced  by  a factor  of 
four  through  the  addition  of  the  0.4%  hydrogen  impurity.  This  result  is  quite 
different  from  the  argon- xenon  mixtures  used  by  Kelly  (Ref.  12) . Kelly  showed 
that  the  addition  of  0.1 % and  0.48%  xenon  to  the  test  gas  argon  did  not  sub- 
stantially change  the  relaxation  length  from  that  for  pure  argon.  The  reason 
is  that  the  excitational  cross-section  for  pure  xenon  is  smaller  than  that  for 
pure  argon.  In  addition,  the  mass  of  atomic  hydrogen  is  markedly  smaller  than 
that  of  the  argon  atom.  These  two  characteristics  of  the  hydrogen  impurity  in 
the  argon  test  gas  give  rise  to  a significant  reduction  of  the  relaxation  length 
or  time. 

3.6  Discussion 

The  following  conclusions  can  be  made  from  the  present  calculations 
and  their  comparison  with  the  experimental  results  for  the  relaxation  processes 
occurring  behind  a strong  shock  wave  moving  in  pure  argon  or  argon  with  a small 
hydrogen  impuri ty: 
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1.  The  initial  process  of  ionization  behind  a strong  shock  moving  in  pure 
argon  is  due  to  atom-atom  collisions.  As  the  number  of  electrons 
increases  sufficiently,  the  process  is  controlled  by  atom-electron 
collisions,  and  subsequently  by  three-body  recombination  (ion- electron- 
electron)  . 

2.  The  new  value  of  the  excitational  cross-section  constant  for  argon-argon 
atom  collisions,  S|r_^r  = 3*5  x 10" 20  cm2/eV,  as  determined  from  a 
comparison  of  predicted  and  measured  data,  is  smaller  than  the  value  of 
1.2  x 10  1 cnr/eV  determined  by  Kelly  (Ref.  10)  and  slightly  larger 
than  the  value  2.5  x 10"2°  cn^/eV  given  by  McLaren  and  Hobson  (Ref.  l4) . 

3.  Owing  to  the  low- mass  of  the  hydrogen  atom  and  the  large  excitational 
cross-sectioh  between  hydrogen  atom-atom  and  hydrogen  atom-electron 
collisions,  the  itpurity  effect  of  hydrogen  in  the  argon  test  gas  is 
to  markedly  reduce  the  relaxation  length  or  time  of  the  ionization 
processes  behind  the  shock  wave,  even  when  the  impurity  level  of 
hydrogen  is  as  lov;  us  0.4$. 

4.  The  effects  of  radiation  losses  and  of  the  wall  boundary  layer  are 
problems  for  further  study. 


4.  LAMINAR  BOUNDARY-LAYER  FLOW  OF  PARTIALLY  IONIZED  ARGON  BEHIND  A STRONG 
SHOCK  WAVE 


4.1  Introduction 

The  nonlinear  partial  differential  equations  for  most  boundary- layer 
problems  are  difficult  to  solve.  As  a result,  many  researchers  have  resorted 
to  using  simplifying  similarity  transformations.  In  the  similarity  approach 
the  system  of  partial  differential  equations  is  reduced  to  a system  of  ordinary 
differential  equations,  which  can  usually  be  solved  by  standard  integration 
techniques.  These  so-called  similarity  solutions  are,  however,  limited  to 
certain  types  of  flows.  For  this  reason,  and  because  of  the  mathematical 
difficulties  encountered  in  obtaining  exact  solutions  for  general  boundary 
layer  cases,  approximate  method"  have  also  been  developed. 

One  approximate  method  of  interest  for  this  work  is  based  on  Harman' s 
momentum  integral,  which  ha,s  been  ext  ided  by  Pohlhausen.  This  approach  is 
now  known  as  the  Karman-Poh: hausen  integral  method.  For  such  integral  methods 
certain  assumptions  are  made  as  to  the  form  of  the  unknown  functions,  which 
reduces  the  problem  to  solving  a set  of  ordinary  differential  equations.  By 
satisfying  appropriate  boundary  conditions  at  the  wall,  the  velocity  prr^ile 
normal  to  the  wall  through  the  boundary  layer  is  reduced  to  a function  of 
only  one  independent  variable.  Libby  and  Morduchow  (Refs.  21  and  22)  have 
extended  the  Karman-Pohlhausen  method  to  a sixth-degree  velocity  profile  and 
a seventh-degree  s tagnation- enthalpy  profile . Chung  and  Anderson  (Refs.  23, 

24)  have  successfully  applied  this  method  to  the  boundary  layer  flow  of  a 
dissociated  gas  over  a flat  plate.  Another  method  developed  by  Dorodnitsyn 
(Ref.  25)  has  been  applied  by  Pallone  (Ref.  26)  to  solve  a mass- transfer 
problem,  and  also  applied  by  Lo  (Ref.  27)  to  solve  a corner-flow  problem. 

This  latter  method  combines  the  Dorodnitsyn  integration  scheme  and  Karman- 
Pohlhausen  integral  method.  Note  also  that  this  approach  was  used  by  Pallone, 
Moore  and  Erdos  (Ref.  28)  in  solving  the  boundary  layer  equations  for  dissoci- 
tion  and  ionizing  air  in  a nonequilibrium  flow. 
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Another  powerful  method  of  solving  the  boundary  layer  equations  is 
the  finite-difference  approach.  Blottner  (Ref,  29)  has  applied  this  numerical 
method  in  solving  nonequilibrium  laminar  boundary  layer  flows  of  an  ionized 
gas.  Kendall  and  Bartlett  (Ref.  30)  have  combined  the  finite- difference 
approach  and  matrix- inversion  techniques  (so-called  integral-matrix  method) 
to  obtain  a nonsimilar  solution  of  a multicomponent  laminar  boundary  layer. 

The  laminar  boundary  layer  behind  a shock  wave  moving  in  a perfect 
gas  has  been  studied  quite  extensively  (Ref's.  31  to  3*0 . The  correct  boundary 
layer  equations  were  used  in  Ref.  32.  Values  of  skin  friction  and  heat  trans- 
fer coefficients  were  obtained.  However,  velocity  and  temperature  profiles 
through  the  boundary  layer  were  not  given.  Mirels  (Refs.  33  and  3*0  solved 
the  laminar  and  turbulent  boundary  layer  equations  for  the  flow  behind  a 
shock  wave.  The  most  extensive  calculations  for  laminar  boundary  layers 
behind  a strong  shock  wave  for  a dissociating  gas  cai  be  found  in  Refs.  35 
to  37.  Khoos  (Ref.  38)  studied  the  Rayleigh  and  shock- tube  end-wall  boundary 
layers  for  partially  ionized  argon.  In  his  analysis,  the  basic  assumption 
is  that  of  thermochemical  equilibrium.  Fay  and  Kemp  (Ref.  39)  considered 
the  heat  transfer  to  a shock- tube  end  wall  from  an  ionized  monatomic  gas  for 
both  frozen  and  equilibrium  flows; 

For  the  steady  laminar  boundary  layer  in  a partially  ionized  gas, 
many  researchers  (Refs.  4o  to  42)  have  studied  the  effects  of  a nonequilibrium 
electron  temperature.  Back  (Ref.  43)  solved  this  problem  based  on  a similarity- 
solution  approach  for  a frozen  low- speed  flow.  The  effects  of  high  temperature 
on  transport  properties  and  nonequiiibrium  flow  remain  unsolved. 

For  the  present  theoretical  work  the  boundary  layer  is  assumed  to  be 
laminar.  Note  also  that  the  boundary  layers  in  the  actual  experiments  (Refs. 

15  and  18)  are  assumed  to  be  laminar. 

4-2  Transport  Properties  of  Partially-Ionized  Argon 

The  kinetic  theory  of  gases  provides  a means  of  estimating  the 
transport  coefficients  of  a partially-ionized  gas.  In  this  section,  transport 
properties  of  parti ally- ionized  argon  gas  are  considered,  as  based  on  the 
mixture  rule  of  Fay  and  Kenp  (Ref.  39). 

For  the  mixture  of  atoms,  ions  and  electrons  the  thermal  conductivity 
(k)  can  be  calculated  by  means  of  the  approximate  mixture  rule; 
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The  respective  symbols  kj  and  Xj  are  the  thermal  conductivity  and  mole  fraction 
of  the  pure  component  j,  and  is  the  effective  hard-sphere  cross-section  for 
a collision  pair  j and’i.. 


The  thermal  conductivities  of  argon  atoms  ( k^,)  and  of  a completely 
singly-ionized  gas  (ks)  are  given  by  the  following  approximate  expressions: 


KAr  = 5,8  x -*-0"7  t3//1+ 


cal/cm  sec°K 


= 4.4  x 10"13  T.  /in  A cal/cm  sec°K 


(4.3) 
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The  ratio  of  the  Debye- distance  to  the  impact  parameter  (for  90°  deflection)  is 
denoted  by  A,  which:  <j  given  by  the  following  relation: 

A = 1.24  x 104  I3/2/n  1/2 
. ' e 

where, T is  in  °K  and  n<a  in  cm-3. 

The  thermal  conductivities  for  the  electron ’gas  (k6)  and  argon  ions 
(KAr+)  can  be  related  to  (ks)  as  shown  below: 


k = (1  + */2)  K 
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From  Eqs.  4.1  to  4.6  the  thermal  conductivity  (k)  for  partially  ionized  argon 
is  given  below: 
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The  viscosity  of  partially  ionized  argon  can  be  calculated  from  simple 
kinetic  theory,  and  it  is  given  below: 
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where 


u = r 8-~—  1 

L ^ V j 


is  the  mean  thermal  speed  of  the  argon  atoms.  The  electrons  make  no  contribution 
to  the  viscosity  because  of  their  extremely  low  mass. 
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The  ambipolar  diffusion  coefficient  Dg^  is  defined  in  terms  of  the 
atom-ion  diffusion  coefficient  D^, 
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(4.9) 


From  kinetic  theory  Da^  is  related  to  the  ratio  of  the  drive  velocity  and  the 
electric  field  E (Ref.  30) , 
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The  ion  temperature  Ty  may  be  taken  as  that  for  which  the  mean  thermal  speed 

r 8kT  -jl/2 

equals  the  drift  velocity,  — . From  experimental  data,  the 

L mAr  J 

effective  hard-sphere  collision  momentum- exchange  cross-sections  for  argon 
atom-atom  collisions  (Qea),  and  for  argon  ion-atom-collisions  (Qai)  can  be 
expressed  as  functions  of  the  gas  temperature  (Refs.  -l8  and  19), 
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Qai  ~ 1.44  T0*16  Qaa 


(4.11) 


The  effective  hard-sphere  argon  ion-ion  collision  cross-section  is  quite 
large  ( ~ 103  to  104  A2);  for  example,  see  Fig.  7 of  Ref.  38. 

Finson  and  Kemp  (Ref.  44)  have  assumed  that  the  pure  ion  viscosity 
can  be  obtained  from  the  ion's  thermal  conductivity  by  using  monatomic  perfect- 
gas  relations,  and  they  developed  an  expression  for  the  Prandtl  number, 


where 


(4.12) 


and 
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1 - a 
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The  degree  of  ionization  for  an  argon  plasma  with  frees bream  and 
surface  conditions  given  in  Table  2 is  plotted  in  Fig.  12  as  a function  of  gas 
temperature  (T)  for  an  equilibrium  flow.  Predictions  for  the  Prandtl  number 
(Pr) , Lewis  number  (Le)  and  the  density-viscosity  product  ratio  (C  = pn/pe!-ie, 
where,  "e"  denotes  the  boundary  layer  edge),  all  obtained  from  the  preceding 
equations,  are  shown  in  Fig.  13 . These  results  for  argon  were  specifically 
based  on  the  initial  conditions  summarized  in  Table  3,  and  the  tenperature 
range  of  300  to  l4,000°K  covers  the  surface  wall  to  frees tream  tenperature 
range  for  the  boundary  layer.  At  the  wall  where  T is  approximately  300°K, 

Pr,  Le  and  C are  equal  to  2/3,  0.44  and  I,  respectively.  At  the  other'  end  of 
the  temperature  range  (l4,000°K  and  higher)  the  gas  is  significantly  ionized, 
and  these  quantities  become  small.  For  a completely  ionized  gas  (T  > 20,000°K) 
the  thermal  conductivity  is  due  almost  entirely  to  electrons,  and  the  electrical 
conductivity  is  strongly  dependent  on  tenperature.  For  these  reasons  the  Prandtl 
and  Lewis  numbers  are  very  small  at  high  temperatures . Note  that  for  a completely 
ionized  gas  or  plasma,  the  viscosity  given  by  Eq.  4.8  may  be  related  to  the  elec- 
tron thermal  conductivity  given  by  Eq.  4.5  (take  the  limit  as  na  -» 0) . Then, 


and  the  Prandtl  number  is  of  the  order  of  l/lOO  for  a completely  ionized  plasma. 
4.3  Karman-Pohlhausen  Integral  Method 

Since  the  original  work  of  Karman  and  Pohlhausen,  their  momentum- 
integral  method  for  solving  boundary- layer  problems  has  been  studied  and 
used  by  many  researchers.  Libby,  Morduchow  and  Bloom  (Ref.  45)  and  Morduchow 
(Ref.  22).  have  made  a critical  study  of  the  integral  method  for  a compressible 
laminar  boundary  layer  for  a perfect  gas.  Chung  and  Anderson  (Ref.  24)  have 
applied  this  method  to  the  boundary  layer  over  a flat  plate  for  a dissociating 
gas.  The  basic  principle  of  the  integral  method  is  to  reduce  the  set  of  non- 
linear partial  differential  equations  for  the  boundary  layer  to  a set  of  ordinary 
differential  equations.  Additionally,  the  profile  through  the  boundary  layer  for 
each  flow  or  dynamic  variable  is  represented  by  a polynomial  equation  of  appro- 
priate degree. 

In  this  section,  the  integral  method  is  formulated  in  a general 
fashion  for  a steady,  two-dimensional  laminar,  compressible,  boundary- layer 
flow  of  a reacting-gas  mixture.  The  freestream  flow  outside  the  boundary 
layer  is  assumed  to  be  in  equilibrium.  It  is  also  assumed  that  the  electron 
and  atom  temperatures  are  equal  in  the  boundary  layer,  or  there  is  one  charac- 
teristic temperature,  T.  If  the  wall  velocity  is  zero,  we  have  the  usual 
boundary  layer  behind  a moving  shock  wave,  we  let  the  wall  velocity  be  equal 
to  the  negative  value  of  the  shock  velocity  ur  attach  the  coordinate  system 
to  the  shock  front,  as  a matter  of  mathematical  convenience.  By  taking  the 
x-axis  along  the  wall  surface  and  the  y-axis  perpendicular  to  the  wall,  the 
boundary- layer  equations  for  a partially-ionized- argon  flow  of  either  type 
are  given  b'y  Eqs.  2.19,  2.20,  2.30  and  2.31. 

The  boundary  conditions  for  the  boundary- layer  flow  behind  a 
moving  shock  wave  are  listed  below, 
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In  the  above  expressions  the  following  conditions  were  used: 
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Equations  4.24  to  4.26  are  the  basic  integral-differential  equation"  for  the 
boundary- layer  flow.  Approximate  solutions  to  these  equations  can  be  obtained 
by  assuming  f , g and  a are  functions  of  only  T)  , which  satisfy  certain  conditions 
at  rj  = 0 and  T)  = 1.  These  boundary  conditions  will  be  discussed  later. 

The  general  solution,  based  on  the  previous  integral-differential 
equations,  can  be  obtained  for  the  special  case  of  a flow  without  a pressure 
gradient.  For  this  special  case  the  simplified  integral-differential  equations 
for  the  flow  over  a flat  plate  are  given  below, 
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The  boundary  conditions  for  velocity,  enthalpy  and  degree  of  ionization 
are  given  below;  • 
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This  last  relation  will  he  developed  subsequently. 
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where,  the  prime  denotes  differentiation  with  respect  to  T).  The  edge  of  the 
sheath  is  now  the  boundary  condition  for  the  boundary  layer.  Consequently,  K, 
the  catalytic  efficiency  coefficient,  is  determined  as  follows:  The  wall  is 
at  a floating  potential.  The  momentum  equation,  Eq.  4.32a,  can  be  obtained 
(Ref.  4l)  from  Langmuir-probe  theory,  if  the  electrons  have  a Maxwellian 
velocity  distribution, 
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where 
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and  &t>  is  the  potential  between  the  wall  and  plasma,  e denotes  the  electronic 
charge  and  s denotes  the  sheath  edge.  The  second  relation  (Eq.  4.32b)  is 
obtained  from  the  continuity  of  mass  flow  of  ions  at  the  outer  edge  of  the 
sheath. 
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The  final  relation  (Eq.  4.32c)  is  obtained  (Refs.  4l  and  42)  from  the  continuity 
of  electron- energy  flux  through  the  outer  edge  of  the  sheath, 
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where,  is  the  thermal  conductivity  of  electrons,  V&e  is  the  drift  velocity 
of  electrons,  and  he  is  the  enthalpy  of  electrons.  * 

As  already  mentioned  in  Section  2.4,  the  sheath  is  very  thin  compared 
with  the  boundary- layer  thickness  for  the  present  problem,  and  the  boundary 
conditions  are  taken  to  be  those  at  the  wall.  From  the  above  equations  for  the 
sheath,  z’(x,0)  can  be  related  to  z(x,0)  and  K are  given  by  (Refs.  4l  and  42), 
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The  dimensionless  velocity,  total  enthalpy  and  the  degree  of  ionization, 
defined  by  Eqs.  4.22,  can  be  approximated  by  the  following  polynomial  equations 
for  their  profiles, 
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Morduchow  (Ref.  22)  has  shown  that  sixth  and  seventh  degree  polynomial  equations 
for  the  velocity  and  enthalpy  profiles,  respectively,  appear  to  be  the  most 
promising  approximation.  Chung  and  Anderson  (Ref.  24)  have  assumed  a fifth- 
degree  polynomial  for  the  profile  of  the  degree  of  dissociation  in  a reacting 
dissociative  gas,  and  found  this  approximation  to  be  reasonable,  as  did  others. 


The  other  supplementary  boundary  conditions  for  a^,  bn  and  cn  can  be 
found  by  estimating  Eqs.  2.20,  2.30  and  2.31  and  their  derivatives  with  respect 
to  coordinate  y.  The  following  supplementary  boundary  conditions  for  a boundary- 
layer  flow  without  a pressure  gradient  and  at  i)  = 0 (or  y = 0)  yield, 
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where,  C = pp/p epe  and  the  prime  denotes  differentiation  with  respect  to  t). 
Equations  4.31  and  4.37  produce  seven,  seven,  and  five  boundary  conditions  for 
a^,  bn  and  cn,  respectively. 
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An  esti  ate  of  the  effect  of  C varying  a ross  the  boundary  layer  was 
studied  by  Back  and  Witte  (Ref.  46)  from  the  stagn«cion-point  heat- transfer 
predictions  for  a high- temperature  un-ionized  gas.  Lees  (Ref.  47)  has  set 
Pw^w  = PePe  transfer  calculation.  This  approximation  might  be 

reasonable  for  an  accelerated  ionized  gas  flow  over  a highly  cooled  wall.  For 
simplicity,  we  adopt  the  method  of  Lees  by  assuming  that  C is  constant. 


From  the  given  boundary  conditions,  an,  bn  and  rn  are  listed  below, 
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b5  = -84(1  - b ) + 45bi  + 20b 2 + 6b3 
b = 70(1  - b ) - 36bi  - 15b2  - 4b3 
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c3  = 10(1  - cQ)  - 6ci  - 3c2 
C4  = -15(1  - cQ)  + 8ci  + 3c2 
c 5 = 6(1  - co)  - 3cx  - c2 


The  remaining  two  unknowns  c0  and  hi  can  be  obtained  by  solving  Eqs.  4.29  and 
4.30. 
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where , 


By  using  the  above  results , the  following  relations  can  be  obtained, 

Fl  = (1  - r)  H9  + §07  r(l  ' r)  (4*39a) 

?2  = B0  - Bi  bi  - B2  b2  - B3  b3  (4.39^) 

F3  = c0  - Cl  Cl  - C2  C2  (4.39c) 
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Substitution  of  Fi  into  Eq.  4.28  yields  the  following  expression  for  A , 
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where  ai  = 2(1  - r) . 

The  solution  for  A can  be  obtained  from  Eq.  4.40  by  using  the  initial 
conditions  A = 0 at  I = 0,  as  shown  below, 

A = A0i  (4.42) 

This  solution  implies  that  the  boundary  layer  thickness  is  proportional  to  the 
square  root  of  I (or  x). 

The  boundary  layer  displacement  thickness  5*,  which  is  defined  by  the 
following  expression, 


6* 


1 - 


Pu 

P u 
e e -l 


dy 


(4.43) 


32 


can  be  related,  to  A as  follows: 


5*  = IlL  j 

R<4 


where,  J is  defined  by  the  egression, 


J\f 


(4.44) 


When  Eqs.  4.38  to  4.4o  are  used  in  conjunction  with  Eqs.  4.29  and 
4.30,  the  following  expressions  for  bi(|)  and  c0(|)  are  obtained, 

+ [ \ Bl  + WT0  j bl  ^ \ B°  " \ B2b2  " I B3bs 

, db2  db3  R , ,N  d + 1 RTwK  „ , 

- B2|  w - B3  w - (Le  - 1)  -li Ci  ( 


(4.45) 


5 (1  + r)  + Cx  K A0|  j + J (1  + r)  + | Ci  K 'TKJ1 

+ ^ K ] c0  = b (1  + r)  ‘ \ C£°2  ■ C2e  tr + cnrrj  £ dT1  (4*46) 

u J e e oJ0 

The  derivatives  of  b2,  b3  and  c2  with  respect  to  I are  found  from  Eqs.  4.38. 

The  initial  conditions  for  Eqs.  4.45  and  4.46  are:  bi  is  finite  and 
cQ  = 1 at  | = 0.  These  conditions  will  be  discussed  later  for  the  equilibrium, 
frozen  and  nonequilibrium  flows. 

The  solutions  for  the  velocity,  total  enthalpy  and  degree  of  ioniza- 
tion profiles  can  be  obtained  from  Eqs.  4.34,  4.35  and  4.38  and  from  the  results 
of  Eqs.  4.45  and  4.46. 

The  use  of  the  integral  method  to  obtain  a solution  for  the  laminar 
boundary  layer  equations  is  widely  accepted  in  aerodynamics..  Generally, 
solutions  are  determined  by  integrating  the  appropriate  equations  in  two 
ways:  first,  in  the  direction  normal  to  the  wall,  using  assumed  polynomial 

distributions;  and,  second,  in  the  flow  direction,  where  new  dependent  vari- 
ables arise  as  a result  of  the  first  integration.  Such  a method  has  proved 
to  be  of  considerable  value  in  assessing  skin- friction  and  heat- transfer 
effects.  However,  such  methods  are  of  limited  usefulness,  since  they  give 
no  indication  of  the  detailed  nature  of  the  flow. 

It  can  be  seen  from  Eqs . 4 .24  to  4 .26  that  the  solutions  are  independ- 
ent of  the  variations  of  the  Prandtl  and  Lewis  numbers  inside  the  boundary 
layer,  and  they  are  dependent  only  on  these  numbers  evaluated  at  the  wall.  As 
a result,  the  successful  method  of  Pallone  (Ref.  26),  who  included  the  Dorod- 
nitsyn  (Ref.  25)  integration  scheme  with  the  Karman-Pohlhausen  approach  was  used. 
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The  method  is  formulated  as  follows.  The  viscous  domain  from  the  body  surface 
to  the  outer  edge  of  the  boundary  layer  is  divided  into  N curvilinear  strips. 
The  boundary  layer  equations  can  then  be  integrated  along  a coordinate  normal 
to  the  body  from  the  surface  to  the  boundary  of  each  strip.  The  essential 
difference  between  the  Karman-Pohlhaus en  and  Dorodnitsyn  integral  methods  is 
that,  in  the  former,  the  partial  differential  equations  are  satisfied  only 
on  the  average  over  the  entire  boundary- layer  thickness  and  in  the  latter  the 
partial  differential  equations  are  satisfied  on  the  average  over  each  strip  of 
the  boundary  layer.  Note  that  for  the  laminar  boundary- layer  flow  of  a non- 
equilibrium ionized  gas , the  use  of  the  Dorodnitsyn  integration  scheme  makes 
the  calculations  very  complicated. 

4.4  Equilibrium  Flow  of  Partially-Ionized  Argon 

The  Rayleigh  and  shock- tube  wall  boundary  layers  for  an  e quilibrium 
flow  of  ionized  argon  were  studied  quite  extensively  by  Khoos  (Ref.  38) • He 
concluded  that,  for  a frees tream  temperature  of  12,000°K  or  greater  and  a 
frees tream  pressure  of  the  order  of  1 atm,  the  equilibrium  assumption  is 
typically  correct  for  the  region  of  the  boundary  layer  where  T > 11,000°K,  but 
breaks  down  for  T < 9j000°K,  at  times  of  the  order  of  10  psec  for  a Rayleigh 
boundary  layer,  and  at  a distance  (x)  of  about  5 cm  for  the  shock-tube  boundary 
layer.  Hence,  near  the  frees tream  edge  of  the  boundary  layer  the  flow  is  in 
equilibrium,  and  near  the  wall  surface  the  flow  is  in  nonequilibrium. 


In  the  case  of  ionization  equilibrium  the  degree  of  ionization  (a) 
is  normally  considered  as  a function  of  temperature  (T)  and  pressure  (p) , 
while  the  pressure  is  related  to  a,  T end  density  p by  the  equation  of  state 
(Eq.  2.23).  The  equation  for  a for  an  equilibrium  boundary  layer  is  given 
below  (Eq.  3«l6),  and  is  called  the  Saha  equation  (Ref.  3)> 


a,  = 


RPITI 


)5/2  e<Tl/T>  + 1 


-1/2 


(4.47) 


Note  that  Pj  = 150.27  gm/cm3  and  Ti  = l82,850°K  for  argon  (Table'  3j  > 


Cp  = | R(1  + a)  + | (1  - a2)  R 


(4.48) 


The  basic  equations  for  an  equilibrium  boundary  layer  flow  of  partially 
ionized  argon  are  given  by  Eqs.  2.19,  2.20  and  2.21,  where  a is  related  to  T 
through  the  Saha  equation  (Eq.  4.47).  Thus,  8a/8y  can  be  related  to  8x/8y. 


8a  8a  8t 


(4.49) 


The  supplementary  boundary  conditions  given  by  Eqs.  4.37c  and  4.37&  are  now 
given  as  follows: 


(Le  - 1)  B g 
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J. 


(4.50a) 
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[ fe  8’  (Le  - 1)  P g'j^  = - ^ [ 0 (l  - ±j  ) (4.50b) 
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6 = | RT/l 


(4.504) 


The  variations  of  € and.  p with  T are  shown  in  Fig.  l4.  For  a wall  temperature 
of  300°K  e and  p are  approximately  zero.  Also,  pw  ~ 0 and  p^-  ~ 0.  If  Pr,  Le 
and  C are  assumed  constant  at  the  wall  surface,  Eq.  4.50a  reduces  to, 


and  Eq.  4.50b  becomes. 


2 ^ (Pr  - 1)  (1  - r)2 
1 + (Le  - 1)  p^ 


b3  ~ 0 


(4.51) 


(4.52) 


below. 


The  differential  equation  for  bi  (Eq.  4.45)  can  be  expressed  as  shown 


IBi  lr + I Bi  + iix"  bi  = \ Bo  ■ i Bshs 
L o j 


(4.53) 


The  initial  conditions  for  bi  is  that  bi  must  be  finite  at  1 = 0.  Equation 
4.53  then  yields  the  following  solution, 


r:  \ B - \ B2b2 

Ul  = i-2 § 

1 , R 

2 Bi  PrF 

o 


(4.54) 


In  Refs'.  22,  23  and  24  it  was  shown  that  in  an  equilibrium  boundary 
layer  a similar  solution  can  exist,  as  the  solution  for  f,  g and  a are 
independent  of  5.  However,  if  C is  not  assumed  constant,  then  a similar 
solution  does  not  exist  in  the  general  case. 

4.5  Frozen  Flow  of  Parti ally- Ionized  Argon 

In  the  case  of  a frozen  flow  of  partially- ionized  argon  it  is  assumed 
that  Cl  remains  constant.  Then,  the  temperature  derivative  of  a (i.e.,  2a/2T)  is 
zero.  The  specific  heat  at  constant  pressure  for  a frozen  flow  is  given  by, 


cp  - I B(1  + 0!) 


(4.55) 


If  the  degree  of  ionization  is  constant,  there  is  no  net  production  of  ions  or 
atoms  and  w is  zero.  Hence,  the  electron  conservation  of  energy  equation  (Eq. 
2.30)  reduces  to, 
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The  degree  of  ionization  for  the  frozen  boundary- layer  flow  can  be 
found  using  Eg..  4.58  and  4.36.  The  initial  condition  required  is  that  c0  = 1 
at  I = 0,  since  the  degree  of  ionization  has  its  freestream  value  at  the  leading 
edge  (I  = 0) . After  the  degree  of  ionization  a has  beerufound,  bi  can  be  ob- 
tained from  Eq.  4.57  by  using  the  initial  condition  that  bx  should  be  finite 

at  £ — 0}  /-T  4.  ,5  pm  N/-y 

1„  1 „ „ db3  B u 2 e 
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(4.60) 


j One  can  see  from  Eq.  4.58  that  the  degree  of  ionization  is  dependent  on  £,  and 
therefore  the  profile  for  the  degree  of  ionization  in  the  frozen  boundary  layer 
is  not  similar,  that  is,  a similar  solution  does  not  exist.  The  profile  of  the 
total  enthalpy  (g)  is  still  a function  of  i and  it  is  also  not  similar,  in 
contrast  to  the  equilibrium  results . However,  the  velocity  distribution  for 
I the  frozen  boundary  layer  is  similar  in  the  present  work.  If  the  value  C was 

t not  assumed  to  be  constant,  then  a nonsimilar  velocity  profile  would  also 

| exist . 


4.6  Nonequi librium  Flow  of  Partially- Ionized  Argon 

In  the  nonequilibrium  boundary  layer  the  mass-production  rate  for 
| electrons  (or  ions)  is  not  zero,  but  it  is  a function  of  temperature  and 

; 
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density.  In  order  to  solve  Eqs . 4.45  and  4.46,  the  term  w/p  must  be  evaluated. 
The  colli sional  ionization  rates  for  a partially-ionized  argon  flow  have  been 
discussed  in  Chapter  3*  The  relationship  between  w/p  and  the  reaction  rate 
dCf/dt  given  by  Eq.  3«l4  is, 


u 

P 


+ a 

e 


(4.61) 


where  da  and  de  are  the  atom- catalyzed  and  electron- catalyzed  reaction  rates, 
respectively,  as  described  in  Chapter  3. 


In  the  present  theory  it  is  assumed  that  the  electron  and  atom  tempera- 
tures are  equal  in  the  boundary  layer,  or  one  temperature  exists.  Consequently, 
da  and  aQ  are  functions  only  of  this  gas  temperature.  The  gas  temperature  in 
the  boundary  layer  is  expected  to  be  less  than  l6,000°K  for  a freestream  tempera- 
ture (T  ) of  13,000°K.  For  this  temperature  range  the  significant  contribution 
to  the  mass -production  rate  w/p  is  due  to  ion-electron-electron  recombination, 
and  the  contributions  due  to  atom-atom  and  atom-electron  collisions  should  be 
small. 


The  initial  conditions  imposed  for  bi  and  c0  at  t = 0 are  the  same  as 
for  the  case  of  a frozen  boundary- layer  flow.  Thus,  bi  is  given  by  Eq.  4.60 
and  c0  = 1. 

4.7  Comparison  of  Theoretical  and  Experimental  Results 
4.7.1  Shock-Tube  Side-Wall  Boundary  Layer 

For  the  case  of  a laminar  boundary  layer  induced  behind  a moving  shock 
wave,  the  boundary  conditions  for  the  freestream  flow  and  the  wall  are  listed 
in  Table  4.  These  initial  conditions  used  for  the  boundary  layer  calculations 
correspond  to  the  experimental  conditions  for  Brimelow's  experiments  (Ref.  15), 
facilitating  the  comparison  of  theoretical  and  experimental  results.  Note  that 
appropriate  constants  for  argon  used  in  the  calculations  are  listed  in  Table  4. 

The  predicted  equilibrium,  frozen  and  nonequilibrium  density  (Pq/ p)  and 
degree  of  ionization  (a)  profiles  for  an  assumed  laminar  boundary  layer  behind 
a shock  wave  are  shown  in  Fig.  15.  The  initial  conditions  are  for  Case  I, 

Table  4.  The  experimental,  results  from  Ref . 15  for  3*75  cm  behind  the  wave  are 
also  shown.  The  three  predicted  density  profiles  are  very  similar,  and  they  are 
in  good  agreement  with  the  experimental  data.  For  the  three  predicted  degree 
of  ionization  profiles,  the  equilibrium  profile  differs  markedly  from  the  very 
similar  frozen  and  nonequilibrium  profiles,  which  are  in  good  agreement  with  the 
experimental  data. 

For  the  equilibrium  boundary  layer  the  degree  of  ionization  attains 
its  equilibrium  value  instantaneously  at  the  local  temperature  and  pressure. 
Since  the  wall  temperature  is  very  low  a equals  zero  at  the  wall  (y  = 0) . 
However,  it  can  be  shown  by  using  Eq.  4.47  that  (23/c>y)w  is  almost  zero,  or 
(Sy/Sa)w  is  almost  infinite.  Consequently,  the  equilibrium  solution  for  the 
degree  of  ionization  differs  markedly  near  the  wall  from  the  frozen  and  non- 
equilibrium solutions.  Note  that  at  the  freestream  edge  of  the  boundary  layer 
the  equilibrium  solution  for  the  degree  of  ionization  is  in  good  agreement  with 
the  experimental  results,  as  expected.  This  behaviour  is  in  agreement  with  the 
conclusions  of  Knoos  (Ref . 38) • 
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For  bhe  frozen  solution,  the  .ionization  reaction  does  not  occur  and 
the  degree  of  ionization  is  controlled  hy  diffusion.  It  is  also  affected  by 
the  surface  recombination  catalycity,  as  described  in  Section  4.3.  From  the 
initial  conditions  for  cQ,  the  degree  of  ionization  profile  at  x = 0 should 
be  constant,  as  a(y)  = %.  For  x greater  than  zero  the  profile  changes  due  to 
diffusion.  As  x increases  from  zero  the  value  of  Q 'v  decreases  from  0.w  and 
approaches  zero  as  x tends  to  infinity.  The  degree  of  ionization  at  the 
edge  of  the  sheath  layer  is  not  zero.  It  can  be  determined  by  using  Eqs. 

4.33  and  4.31a  for  z'(0).  In  Fig.  15,  the  degree  of  ionization  at  the  edge 
of  the  thin  sheath  layer  is  seen  to  be  very  small,  but  it  is  nonzero. 

The  nonequi librium  solution  for  the  degree  of  ionization  is  very 
similar  to  the  frozen  solution.  The  reason  for  this  behaviour  is  that  the 
gas  temperature  within  the  boundary  layer  is  below  15,000°K  and  the  mass 
production  rate  of  electrons  therefore  has  a minor  effect  on  the  degree  of 
ionization* 

For  Case  II  of  Tables  2 and  4,  the  predicted  equilibrium,  frozen  and 
nonequilibrium  solutions  for  the  density  and  degree  of  ionization  profiles  are 
shown  in  Figs.  l6  and  17,  along  with  the  experimental  results.  The  same 
comments  made  for  the  previous  case  apply  to  these  results.  Additionally,  by 
comparing  theoretical  and  experimental  results  it  can  be  seen  that  near  the 
edge  of  bhe  boundary  layer  the  flow  is  in  quasi- equilibrium  due  to  radiation 
losses  that  become  increasingly  important  with  distance  (20  cm).  Nonequilib- 
rium effects  are  important  near  the  wall  region. 

For  Case  III  of  Tables  2 and  4,  bhe  equilibrium,  frozen  and  nonequi- 
librium solutions  for  the  density  and  degree  of  ionization  profiles  are  shown 
in  Figs.  18  and  19.  For  this  weaker  shock  wave  and  lower  ionization  case, 
the  predicted  and  measured  results  are  in  good  agreement.  At  12  cm  and 
lower  temperatures  radiation  does  not  appear  to  be  too  significant. 

For  the  case  of  a steady  boundary  layer  in  a partially-ionized-gas 
flow  over  a flat  plate,  the  integral  method  described  in  Section  4.3  can  be 
applied  by  setting  r equal  to  zero.  The  significant  difference  between  this 
quasi- steady  boundary  layer  and  the  boundary  layer  behind  a shock  wave  is  in 
the  velocity  profile.  In  the  latter  case  the  velocity  of  the  wall  equals 
zero. 

4.7.2  Quasi-Steady  Flat-Plate  Boundary  Layer 

The  freesbream  and  wall  conditions  assumed  for  the  steady  lamih&rc  . 
boundary  layer  over  a flat  plate  in  the  UTIAS  hypersonic  shock  tube  are 
summarized  in  Table  3. 

For  the  case  of  an  equilibrium  boundary  layer.,  bhe  predicted  degree 
of  ionization  and  temperature  profiles  as  a function  of  T|  (t)B  is  defined  in 
Section  4.8)  for  two  particular  Prandbl  numbers  of  1 and  0.7  are  shown  in 
Fig.  20.  It  can  be  readily  seen  that  the  effect  of  bhe  Prandtl  number  (in 
the  energy  equation)  on  the  degree  of  ionization  and  temperature  is  significant. 

Equilibrium  and  frozen  solutions  for  the  temperature  and  degree  of 
ionization  profiles  for  a flat-plate  boundary  layer  are  conpared  in  Fig.  21, 
and  the  associated  density  profiles  are  conpared  in  Fig.  22.  For  both 
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solutions  the  Prandtl  and  Lewis  numbers  were  set  equal  to  unity.  The  equilibrium 
and  frozen  solutions  exhibit  sizeable  differences.  Note  that  the  frozen  solutions 
displayed  in  both  Fig.  21  and  22  were  calculated  for  the  case  when  x approaches 
infinity,  or  when  = 0.  Then,  an  exact  numerical  solution  can  be  obtained 
easily  by  using  the  similarity  assumption. 

For  a flat-plate  boundary  layer  having  the  boundary  conditions  given 
in  Table  3,  equilibrium,  frozen  and  nonequilibrium  solutions  for  the  density 
profile  are  compared  in  Fig.  23a,  and  those  for  the  degree  of  ionization 
profile  are  given  in  Fig.  23b.  For  these  calculations  the  Prandtl  number, 

Lewis  number  and  the  density- viscosity  product,  C,  were  taien  as  0.67,  0.44 
and  1.6,  respectively.  The  three  predicted  profiles  for  the  density  are  all 
similar,  but  their  agreement  with  Whizten's  experimental  data  (Ref.  18)  is  not 
good.  Note  that  the  minimum  value  of  the  density  is  at  the  edge  of  the  boundary 
layer  for  the  predicted  profile,  while  it  occurs  at  y equal  to  0.7  mm  in  the 
experimental  profile.  This  disagreement  is  presently  being  investigated. 


The  predicted  profiles  for  the  degree  of  ionization  are  quite  different 
(Fig.  23b).  Also,  a signiiicant  discrepancy  exists  between  the  theoretical  and 
experimental  results.  The  overshoot  of  the  degree  of  ionization  beyond  its 
freestream  value  is  not  understood  at  present,  Some  of  the  disagreement  between 
the  predicted  and  experimental  results  could  be  due  to  errors  in  the  integral 
method.  As  described  in  Section  4.3,  the  present  Kaiman-Pohlhausen  integral 
method  uses  flow  quantities  that  are  averaged  over  the  entire  boundary  layer 
thickness . 


4.8  Comparison  of  Integral  Method  and  Exact  Numerical  Calculations 

Similar  solutions  and  exact  numerical  calculations  were  made  for 
identical  boundary- layer  problems  for  comparison  purposes.  By  an  exact 
solution  it  is  meant  that  the  Blasius-type  solution  is  obtained  by  using  the 
Howarth  transformation  (Ref.  48)  specified  below. 
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The  subscript  B refers  to  the  Blasius  type  of  variable.  In  order  to  distinguish 
between  the  notation  of  the  dimensionless  velocity  in  this  section  and  in  Section 
4.3,  we  define 
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For  a steady  laminar  boundary  layer  of  a partially-ionized  gas,  Back 
(Ref.  47)  has  discussed  the  similarity  solution  approach  for  a low-speed  flow. 
The  basic  equations  for  the  boundary  layer  (Eqs.  2.19,  2.20,  2.30  and  2.31) 
can  be  transformed  to  yield  the  following  expressions. 
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The  prime  denotes  differentiation  with  respect  to  t]g. 

The  flow  over  a flat  plate  is  considered,  with  Pr,,  Le  and  C held  constant 
for  simplicity.  For  the  equilibrium  case  the  ion  concentration  gradient  is  ex- 
pressed in  terms  of  the  enthalpy  gradient.  The  basic  equations  (Eqs.  4.63  and 
4.64)  reduce  to, 
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The  symbol  0 is  defined  by  Eq.  4.50c.  The  boundary  conditions  for  the  equili- 
brium, similar  boundary- layer  equations  are  as  follows: 

fB(0)  = fg(o)  = 0 
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For  the  case  of  a frozen  boundary- le.yer  flow,  Eqs-.  4.63  to  4.65  take 
the  following  forms, 
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The  boundary  condi bions  are  as  given  in  Eqs.  4.68  with  the  following  additional 


ones , 
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In  order  bo  compare  results  of  both  the  integral  method  and  the  exact 
similar  solutions,  a common  coordinate  must  be  used.  From  Eq.  4.21  we  have 
the  following  expression: 
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Define  tj*  as  the  common  coordinate, 
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(4.74) 


On  the  other  hand,  T|*  can  be  related  to  by, 


n* 


V dnB 


(4.75) 


From  Eqs.  4.74  and  4.75  the  relations  between  the  derivatives  re. 
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o 

where,  0 represents  an  arbitrary  parameter. 

For  a steady  boundary  layer  over  a flat  plate  with  the  frees bream 
and  wall  conditions  as  given  in  Table  2,  the  exact  numerical  profiles  of 
velocity,  enthalpy,  degree  of  ionization  and  temperature  are  displayed  in 
Fig.  24.  The  corresponding  density  profile  is  shown  in  Fig.  25.  Also,  the 
effects  of  a different  Prandtl  number  on  the  density  are  illustrated  in  this 
figure.  The  exact  numerical  temperature  profile  has  also  been  plotted  in 
Fig.  26,  where  it  is  compared  with  the  integral  method  profile.  The  results 
obtained  by  the  two  different  methods  are  in  excellent  agreement. 

For  an  equilibrium  boundary- layer  flow  over  a flat  plate,  exact  numeri- 
cal and  integral  method  results  for  the  degree  of  ionizc.tion  and  temperature 
profiles  were  computed.  These  results  are  shown  in  Fig.  27,  and  it  can  be 
seen  that  both  sets  are  in  excellent  agreement. 
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The  comparison  of  the  exact  numerical  and  integral  method  results 
serves  the  particular  purpose  m showing  the  accuracy  of  using  polynomial 
expressions  in  the  integral  method.  From  the  comparisons  shown  in  Figs.  26 
and  27,  it  can  be  concluded  that  the  present  integral  method  is  in  good 
agreement  with  the  exact  numerical  one. 

4.9  Discussion 

The  integral  method  has  been  used  to  study  both  the  laminar  wall 
boundary  layer  induced  by  a moving  shock  wave  and  the  quasi- steady  flat-plate 
laminar  boundary  t.ayer.  The  following  concluding  remarks  can  be  made. 

1.  The  integral  method  is  in  good  agreement  with  experimental  data  for 
the  laminar  wall  boundary  layer  behind  a shock  wave.  Near  the  edge  of  the 
boundary  Layer  the  flow  is  expected  to  be  in  equilibrium,  and  near  the  wall, 
the  flow  is  expected  to  be  in  a nonequilibrium  or  in  a frozen  state. 

2.  The  fifth,  sixth  and  seventh  order  polynomial  equations  for  the 
degree  of  ionization,  velocity  and  enthalpy  profiles,  respectively,  are 
sufficiently  accurate  for  the  present  method. 

3.  A comparison  between  some  predicted  and  experimental  results  for  the 
quasi-steady  flat-plate  boundary  layer  indicates  that  there  are  discrepancies. 
Part  of  the  disagreement  may  be  due  to  the  fact  that  actual  variation?  of 
Prandtl  number,  Lewis  number  and  the  density  viscosity  product  were  n 
included  in  the  analysis.  A more  satisfactory  method  might  be  to  combine 

the  Karman-Pohlhausen  integral  method  and  „ne  Dorodnitsyn  integration  scheme 
(Ref.  25). 


It  is  possible  that  the  analysis  and  experiment  may  agree  for  the 
side-wall  boundary  layer  where  radiation  effects  close  to  the  shock  front 
are  not  yet  significant.  However,  in  the  case  of  the  flat  plate  quasi-steady 
boundary  layer,  where  measurements  are  made  far  from  the  shock  front  (20-30  cm)., 
radiation  losses  are  inportant.  For  example,  at  shock  Mach  numbers  of  about 
17,  the  frou-ot ream  ionization  has  dropped  from  16 $ to  %.  It  will  therefore 
be  necessary,  as  a first  step,  to  include  the  new  initial  conditions  in  the 
numerical  analysis.  As  a second  step,  it  appears  that  the  boundary  layer 
equations  wou]-.  have  to  be  solved  by  including  the  effects  of  the  radiation 
terms.  Ultirntely,  a finite-difference  scheme  would  have  to  be  developed 
that  would  permit  the  inclusion  of  variable  Prandtl  and  Lewis  numbers  through 
the  boundary  layer.  Perhaps,  when  these  major  modifications  are  made  in  the 
analysis,  better  agreement  with  the  interferometric  results  might  be  expected. 

Nevertheless,  when  one  looks  at  the  results  obtained  by  Brimelow 
(Ref.  15)  for  the  side-vail  boundary  layer  it  is  clear  that  the  agreement  of 
experiment  with  analysis  is  as  good  close  (3*75  cm)  to  the  shock  front  Mg  ~ 17) 
as  for  large  distances  (20  cm).  This  points  to  some  basic  differences  between 
tne  two  sets  of  data  and  the  analyses.  This  problem  is  being  investigated. 


INTERACTION  OF  A LAMINAR  BOUNDARY-  LAYER  FLOW  AND  A CORNER- EXPANSION  WAVE 
IN  P^RTIALLl -IONIZED  ARGON 


5.1  Introduction 


For  an  ideal  and  inviscid  flow  at  supersonic  speed,  the  flow  around 
a sharp  corner  takes  place  through  a stationary  expansion  wave.  Such  a flow 
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is  completely  described  by  the  well-known  Prandtl- Meyer  relations  (Ref.  49). 

The  steady,  inviscid  flow  of  a partially-ionized  gas  around  a sharp,  convex 
corner  is  much  more  conplex  to  describe  analytically,  and  it  has  been  studied 
in  detail  by  Glass  and  Takano  (Ref.  3) • For  a real  gas  flow  a boundary  layer 
is  produced  on  the  wall  surface.  Its  presence  can  lead  to  a considerable 
modification  of  the  results  predicted  by  inviscid  flow  theory  near  the  surface. 
The  interaction  of  a laminar  supersonic  or  hypersonic  boundary  layer  with  a 
steady  corner- expansion  wave  for  a partially-ionized  gas  is  an  important 
problem  for  aerospace  control  surfaces  and  entry  into  planetary  atmospheres. 
This  problem  is  still  not  well  understood  both  theoretically  and  experiment- 
ally. 

The  interaction  of  a laminar  boundary  layer  and  a corner- expans ion 
wave  for  a perfect  gas  has  been  investigated  extensively  (Refs.  27,  50  and  51). 
A comprehensive  solution  to  this  problem  has  been  given  by  Lo  (Ref.  27)  who 
used  Dorodnitsyn's  integral  method  to  solve  the  Na.vier-Stokes  equations.  To 
solve  the  same  problem  Sullivan  (Ref.  51)  suggested  a simpler  approach  called 
the  cold-wall  similarity  method.  This  method  predicts  successfully  the  major 
features  of  the  flow  at  the  edge  of  the  boundary  layer.  The  main  advantages 
of  Sullivan's  method  are  that  the  calculations  are  simple  to  perform  and  the 
results  are  in  good  agreement  with  Lo's.  for  the  pressure  distribucion  and 
boundary- layer  thickness.  Additionally,  the  results  of  Sullivan's  analysis 
are  in  good  agreement  with  experimental  data  (Ref.  52). 

For  an  inviscid  flow,  Glass  and  Takano  (Ref.  3)  made  a detailed 
study  of  nonequilibrium,  frozen  and  equilibrium  expansion  flows  of  ionized 
argon  around  a sharp  comer.  The  numerical  calculation  procedure  was  based 
on  the  method  of  characteristics,  and  it  is  similar  to  that  for  dissociating 
oxygen  (Ref.  53).  Furthermore,  Glass  and  Igra  (Refs.  54,  55)  showed  that 
the  experimental  results  substantiate  the  analysis  of  Gj.ass  and  Takano  for 
a nonequilibrium  supersonic  corner- expansion  flow  of  partially  ionized  argon. 

In  this  chapter,  Sullivan's  method  is  applied  to  solve  the  inter- 
action of  a laminar  boundary  layer  and  a corner- expansion  wave  for  a partially- 
ionized  gas.  For  simplicity  the  assumption  of  a frozen  flow  in  the  boundary 
layer  is  made  from  the  onset.  A complete  solution  for  the  corner- expansion 
problem  would  be  based  on  the  coupled  solutions  for  inviscid  flow  (Ref.  3) 
and  the  Navier-Stokes  equations,  which  makes  the  analysis  complex, 

5.2  Cold-Wall  Similarity  Method  for  a Supersonic  Ionized  Argon  Flow 

The  assumptions  of  hypersonic  flow  and  a simple  expansion  wave  used 
in  Sullivan's  analysis  (Ref.  51)  are  removed,  in  order  to  apply  the  cold-wall 
similarity  method  to  the  case  of  a supersonic  flow. 

The  boundary  layer  displacement  thickness  5*,  defined  by  Eq.  4.43, 
can  be  expressed  as, 
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(5.2b) 


The  coordinate  t)b  is  given  by  Eq.  4.62b.  A function  G(x)  is  defined  as, 


G(x)  = 
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(5.3) 
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The  basic  relationship  between  5*  and  pe/Poo  is, 
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Equation  5.1  can  be  written  in  an  alternate  way  in  terms  of  the  freestream 

1/2 


variables  M^,  p^  and 
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Note  that  z is  proportional  to  x,  and  ‘"hat  7*  equals  5/3  for  the  frozen  flow. 

From  the  cold-wall  similarity  assumption,  G is  independent  of  x,  but 
it  is  of  constant  value.  Differentiating  Eq.  5.5  with  respect  to  z,  then, 
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P5‘  = dP/dZ 

p = Pe/P* 
R =J'p<3Z 


For  a frozen  flow  ae  is  constant  and  the  following  relations  are 
applicable, 
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The  rate  of  growth  of  the  boundary  layer  displacement  thickness 
immediately  downstream  of  the  corner  is, 


* dB* 

dB* 

V 

cbc 

Note  that  aw  is  the  corner  turning  angle.  The  rate  of  growth  of  the  boundary- 
layer  displacement  thickness  downstream  of  the  corner  can  be  related  to  the 
deflection  angle  v(p). 

ir  * “t  • v<p>  <5,i: 

v(p)  = v(Me)  - v(Meo)  (5. 11 

The  Prandtl-Meyer  function  for  a frozen  flow  v(m)  is  given  below. 


v(M)  . 
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Substituting  Eq.  5.13  into  Eq.  5.6  leads  to  the  following  equations  for  the 
boundary- layer  and  exp an si on- wave  interaction  problem, 
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dz  2R 


^R  M ' 
00 
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Solutions  to  these  two  ordinary  differential  equations  using  a Runge-Kutta 
technique  gives  the  pressure  distribution  pe/Poo«  The  initial  conditions  for 
P and  R are  obtained  from  the  values  immediately  upstream  of  the  corner 
(x  = xc);  that  is,  P = 1 and  R = Rc  (see  Fig.  28).  Note  that  the  equation  for 
the  skin  friction  is, 
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(5.18) 


5.3  Model  Calculations 
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The  frees tream  and  wall  conditions  for  the  following  calculations  of 
a laminar  boundary  layer  interacting  with  an  expansion  wave  are  given  in 
Table  2.  Note  that  xc  (Fig.  28)  has  been  set  at  13.5  cm.  The  value  of  G 
(Eq.  5.3)  is  obtained  by  solving  Eqs.  4. 69  to  4.71  (see  Appendix  A),  giving 
G = 0.274(12). 

Calculated  results  of  pressure,  temperature  and  density  distribution 
at  the  edge  of  the  boundary  layer  as  a function  of  distance  are  shown  in 
Figs.  29  to  31j  respectively.  In  each  case  the  results  for  four  different 
turning  angles  of  5,  10,  15  and  20  degrees  are  presented.  Similar  results 
for  Mach  number,  boundary  layer,  displacement  thickness  and  skin  friction 
are  given  in  Figs.  32,  33  and  34,  respectively.  It  can  be  seen  that  as  the 
turning  angle  for  the  corner  increases,  the  surface  pressure,  temperature, 
density,  and  skin  friction  decrease,  whereas  the  displacement  thickness  and 
Mach  number  increase. 

Some  experimental  data  on  the  interaction  of  a boundary  layer  and 
a corner-expansion  wave  are  given  in  Ref.  18.  Profiles  for  the  degree  of 
ionization  and  density  have  been  measured.  As  in  the  flat-plate  quasi-steady 
boundary  layer  (Figs.  23a,  b.) , the  post-corner  boundary- layer  profiles  of  these 
quantities  are  in  disagreement  with  the  present  analysis.  In  this  case  as 
well,  the  reasons  are  being  investigated. 

5.4  Discussion 

The  cold-wall  similarity  method  can  be  used  to  predict  the  flow 
quantities  at  the  boundary  layer  edge  downstream  of  the  corner.  It  is 
believed  thrt  these  predictions  are  quite  reliable.  However,  it  is  doubtful 
that  the  method  can  predict  accurately  the  heat- transfer  and  skin- friction 
coefficients,  owing  to  the  assumption  of  similarity.  For  the  cold-wall 
similarity  method  f"(0)  is  assumed  constant  downstream  of  the  corner.  Actually, 
f"(0)  is  a function  of  d(pe/pw)/dx.  Near  the  corner  this  distance  derivative 
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of  pressure  is  a maximum  and  decreases  to  zero  only  as  x increases  to  infinity. 
Since  the  skin- friction  coefficient  is  proportional  to  f"(0),  the  predicted 
coefficient  is  therefore  only  as  accurate  as  the  underlying  assumption,  which 
is  believed  to  be  inaccurate.  Similarly,  the  heat-transfer  coefficient  would 
not  be  predicted  accurately. 

A new  method,  based  on  an  iterative  procedure,  is  suggested  below. 

1.  Calculate  Pe/pro  and  d(pe/pco)/dx  by  the  cold-wall  similarity  method. 

2.  Calculate  dug/dx  from  the  results  of  step  1. 

3.  Substitute  due/dx  into  Eqs.  4.63  to  4.65  and  solve  them. 

4.  From  the  new  results  for  f ' , g and  g'  calculate  G(x) . 

5.  By  using  the  cold-wall  similarity  method  recalculate  Pe/pM  and  d(pe/pM)/dx. 

6.  Repeat  steps  2,  3S  4,  5 and  6 until  the  resulting  error  is  satisfactorily 
small. 

This  new  method  should  enable  the  flow  quantities  at  the  boundary- layer  edge 
and  on  the  wall  surface  to  be  predicted  as  accurately  as  with  the  finite 
difference  method.  Similarly,  the  velocity,  enthalpy,  and  degree  of  ionization 
profiles  for  the  boundary  layer  downstream  of  the  corner  would  be  predicted 
accurately.  However,  the  advantage  is  that  the  new  method  should  provide  a 
means  of  predicting  accurately  the  heat- transfer  and  skin- friction  coefficients. 

A complete  solution  for  the  flow  field  around  a corner  can  be  obtained 
by  combining  the  boundary- layer  solution  described  in  Chapter  4 with  the  solution 
of  Glass  and  Takano  (Ref.  3)  for  the  expansion  wave.  However,  the  difficulty  of 
defining  a common  coordinate  for  both  parts  of  the  flow  field  - viscous  and 
inviscid  - would  first  have  to  be  overcome. 


6.  CONCLUSIONS 


A study  was  made  of  shock- wave  structure,  laminar  boundary  layers  in 
nonequilibrium  flows  of  partially-ionized  argon  induced  by  a shock  wave  on 
shock  tube  side  walls  and  over  a flat  plate,  and  of  the  interaction  of  a 
laminar  boundary  layer  with  a corner- expansion  wave. 

A new  value  of  the  ionizational  cross-section  constant  for  argon 
atom-atom  collisions,  was  determined  as,  = 3*5  x 10~20  cm2/eV,  by 

comparing  results  from  a theoretical  two-step  collision  mod.:!  with  experimental 
data  of  shock  structure  and  its  relaxation  processes.  It  was  confirmed  that 
the  ionization  rate  is  controlled  by  atom-atom  collisions  initially  and  then 
by  the  more  efficient  electron-atom  collisions.  When  the  electron  number 
density  increases  to  a certain  level,  the  electron-production  rate  due  to 
atom  and  electron  collisions  drops  quickly  and  the  process  reaches  equilibrium. 
Owing  to  the  extremely  low  mass  of  a hydrogen  atom  and  the  large  excitational 
cross-sections  for  hydrogen  atom-atom  and  electron-atom  collisions,  even  a low 
level  of  hydrogen  impurity  can  reduce  the  total  relaxation  length  of  the  shock 
wave  in  argon. 
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The  integral  method  of  Karman  and  Pohlhausen  was  applied  to  study 
the  shock-induced  laminar  wall  boundary  J.tyer  and  the  quasi- steady  flat-plate 
laminar  boundary  layer  in  an  ionized  argon  flow.  Equilibrium,  frozen  and 
nonequilibrium  solutions  were  obtained,  compared  with  each  other,  and  with 
available  experimental  results.  The  experimental  data  show  that  the  flow 
near  the  edge  of  the  shock- induced  boundary  layer  on  a wall  is  in  equilibrium 
but  in  the  wall  region  the  flow  is  in  nonequilibrium.  This  conclusion  is  in 
agreement  with  Khoos'  results. 

For  the  quasi- steady  boundary-layer  flow  over  a flat  plate,  there 
is  disagreement  between  the  analytical  and  experimental  results.  It  is  not 
understood  why  the  analysis  should  agree  with  the  wall  boundary  layer  experi- 
mental data  and  not  with  the  flat-plate  data.  Thxs  discrepancy  may  be  due 
to  the  error  in  using  the  integral  method  for  the  analysis.  The  solutions 
based  on  the  integral  method  are  independent  of  variations  of  the  transport 
properties  with  gas  tenperature  and  degree  of  ionization,  which  is  a serious 
deficiency  of  the  present  analysis.  The  present  method  also  suffers  from 
the  fact  that  the  results  are  averaged  over  the  entire  boundary  layer, 
rather  than  reflecting  rapidly  varying  local  properties.  However,  if  these 
were  the  reasons  then  they  should  be  applicable  to  both  boundary  layers. 

This  is  not  the  case.  Perhaps  the  major  difficulty  arises  from  radiation 
losses  as  discussed  at  the  end  of  Section  4.9. 

The  cold-wall  similarity  method  developed  by  Sullivan  was 
extended  and  applied  to  the  interaction  of  a laminar  boundary  layer  and  a 
corner- expansion  flow  for  partially-ionized  argon.  The  entire  flow  was 
assumed  frozen  both  upstream  and  downstream  of  the  corner.  The  hypersonic 
flow  and  simple  expansion- wave  assumptions  were  removed  for  the  present  case 
of  supersonic  flow.  The  cold-wall  similarity  method  predicts  the  pressure, 
temperature  and  Mach  number  at  the  edge  of  the  boundary  layer,  and  the  dis- 
placement thickness,  but  breaks  down  for  predicting  flow  profiles  inside 
the  boundary  layer.  To  overcome  this  difficulty  a new  method  of  solution 
was  proposed.  Available  post-corner,  boundary- layer  experimental  data  are 
also  in  disagreement  with  the  present  analysis.  The  reasons  why  disagreements 
occur  in  the  case  of  the  quasi-steady  boundary  layers  on  a flat  plate  and 
after  a corner-expansion  interaction  are  presently  under  study,  as  discussed 
at  the  end  of  Section  4.9. 
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TABLE  1 


NUMERICAL  CONSTANTS  USED  IN  THE 


UNIVERSAL  CONSTANTS 

Boltzmann  constant 

Avogadro  constant 

Mass  of  an  electron  (m  ) 

' e 

Partition  function  for  electrons 
CONSTANTS  FOR  ARGON 

Molecular  weight  of  an  argon  atom  (m^) 
Gas  constant  (R^  = k/10^) 

Partition  function  for  atoms 
Partition  function  for  ions 
Ionization  temperature  (Tj) 

First  excitational  temperature  (T*) 
Characteristic  density  of  ionization 

(Pi) 

CONSTANTS  FOR  HYDROGEN 

Molecular1  weight  of  a hydrogen  molecule 
(*iH2) 

Molecular  weight  of  a hydrogen  atom 
(n>H) 

Gas  constant  for  a hydrogen  atom  (Rjj) 
Partition  function  for  atoms 
Partition  functipn  for  ions 
Ionization  temperature  (Tj) 

First  excitational  temperature  (T*) 


CALCULATIONS 


1.3803  x 10'16  erg/°K 

OO 

6.023  x 10  per  mole 

Op. 

9.1066  X 10  g 
2 

39.944  g/mole 
0.20813  x 107  erg/g°K 
1 
6 

l82,850°K 

135,000°K 

150.27  g/cm3 


2.016  g/mole 

1.008  g/mole 
8.3135  x 107  erg/ g°K 
2 
1 

157,000°K 
118, 380 °K 


TABLE  2 

INITIAL  CONDITIONS  FOR  A STRONG  SHOCK  WAVE 
MOVING  IN  A QUIESCENT  GAS 


Case 

Gas 

P^( torr) 

Ms 

T}.(  °K) 

I 

pure  argon 

5.6 

16.5 

298.7 

II 

pure  argon 

5.12 

16.53 

^6.6 

III 

pure  argon 

5-09 

13.59 

296.7 

IV 

argon  + 0.4$  hydrogen 
(by  pressure) 

5.17 

16.68 

297.4 

TABLE  3 

FREESTREAM  AND  WALL  CONDITIONS  FOR  A LAMINAR  BOUNDARY  LAYER 


OVER  A FLAT  PLATE  IN  PURE  ARGON 


Uco 

4778  m/sec 

Mco 

2.4 

2238.4  torr 

*00 

13,393?K 

T 

•w 

298. 2 °K 

M 

s 

16.98 

X 

13.5  era 

b) 


X,  u 
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1 TWO  DIFFERENT  COORDINATE  SYSTEMS  FOR  A BOUNDARY  LAYER  INDUCED  BEHIND 
A MOVING  SHOCK  WAVE,  a)  COORDINATE  SYSTEM  FIXED  TO  THE  WALL,  b) 
COORDINATE  SYSTEM  FIXED  TO  THE  SHOCK-WAVE  FRONT. 
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FIG.  3 SHOCK  STRUCTURE  IN  A NONEQU I L I BR IUM  ARGON  FLOW. 


VARIATION  OF  DEGREE  OF  IONIZATION  (ct)  WITH  DISTANCE  (x)  THROUGH  RELAXATION  REGION 
BEHIND  A STRONG  SHOCK  WAVE  IN  PURE  ARGON  WITH  pi  = 5-6  TORR,  Kg  = 16.5  AND  Ti  = 2Q8.7°K. 


TABLE  2) 


ELECTRON  TEMPERATURE  (Te)  VARIATIONS  WITH  DISTANCE  (x)  THROUGH 
: px  = 5.6  TORR,  Mg  = 16.5  AND  Ta  = 298. 7°K  IN  PURE  ARGON. 


FIG.  8 ELECTRON  NUMBER  DENSITY  (ne)  AND  IONIZATION 
PRODUCTION  RATES  da  AND  ae  WITH  DISTANCE  i'x) 
THROUGH  THE  RELAXATION  REGION  FOR  pi  = 5.6 
TORR,  Mg  = 16.5  AND  Ti  = 298. 7°K  IN  PURE  ARGON. 


FIG.  9 DENSITY  (p)  AND  DEGREE  GF  IONIZATION  (a)  WITH  DISTANCE  FOR  pi  = 5.12  TORR,  Ms  = 16.53  J 
Ti  = 296. 6° K IN  PURE  ARGON. 


CASE  HI  (TABLE  2) 


FIG.  10  DEGREE  OF  IONIZATION  (a)  WITH  DISTANCE  (x)  THROUGH  THE  RELAXATION  REGION  FOR  pi  = 5.09  TORR 
Ms  = 13-59  AND  Tx  = 296. 7°K  IN  PURE  ARGON. 


EQUILIBRIUM  ARGON  PLASMA 
« 2238  torr 


FIG.  12  DEGREE  OF  IONIZATION  (a)  AS  A FUNCTION  OF  TEMPERATURE  (T)  FOR  AN  EQUILIBRIUM  ARGON  PLASMA 


PRANDTL  NUMBER  (Pr) , LEWIS  NUMBER  (Le)  AND  DENSITY-VISCOSITY  PRODUCT  (£)  VERSUS 
TEMPERATURE  (T)  FOR  AN  EQUILIBRIUM  ARGON  PLASMA. 


VARIATION  OF  £ AND  p AS  A FUNCTION  OF  TEMPERATURE  (INITIAL 
CONDITIONS  GIVEN  IN  TABLE  3) . 


DENSITY  RATIO  (pe/p)  VARIATION  THROUGH  A LAMINAR  WALL  BOUNDARY  LAYER  BEHIND  A SHOCK  WAVE 
IN  PURE  ARGON.  THE  INITIAL  CONDITIONS  ARE  GIVEN  IN  FIG.  15,  EXCEPT  x = 20  CM.  THE 
ERROR  BAR  SHOWS  THE  UNCERTAINTY  IN  THE  QUANTITIES  DUE  TO  THE  ERROR  IN  THE  EXACT  WALL 
POSITION  NEAR  THE  WALL. 


CASE  11  (WBLES  a AN0  4) 
“ EQUILIBRIUM 
frozen 

noneouiubr/um 

• EXPERIMENT  (R£e  ,5) 


DENSITY  j^atto  t / / £ 

s r*® 

X “ 20  CM.  JJV-iTIAL  COaiTlows 


FIG.  20  TEMPERATURE  (T/Te)  AND  DEGREE  OF  IONIZATION  (o/oe)  VARIATIONS  THROUGH 
A QUASI-STEADY  FLAT-PLATE  BOUNDARY  LAYER. 


FIG.  22  DENSITY  VARIATION  (p;  THROUGH  A QUASI-STEADY  FLAT-PLATE  BOUNDARY  LAYER. 


2.5 


FIG.  23a  PREDICTED  AND' MEASURED  DENSITY  VARIATIONS  THROUGH  A QUASI-STEADY  FLAT-PLATE  BOUNDARY 


o 02  0.4  0.6  0.8  1.0 


U/Ue,  H/He,  o C/£e  and  T/Te 


FIG.  24  PREDICTED  PROFILES  OF  FLOW  VELOCITY,  TOTAL  ENTHALPY,  DEGREE  OF  5 

IONIZATION  AND  TEMPERATURE  FOR  A QUASI-STEADY  FLAT-PLATE  BOUNDARY  1 
LAYER. 
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FIG.  28  FLOW  PATTERN  OF  THE  BOUNDARY  LAYER  INTERACTION  WITH  A CORNER 
EXPANSION  WAVE. 
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ATURE  VARIATION  WITH  DISTANCE  FOR  A 
RY  LAYER-CORNER  EXPANSION  FLOW. 
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APPENDIX  A:  METHOD  OF  SOLVING  BOUNDARY-LAYER  EQUATIONS 

The  partial  differential  equations  describing  the  boundary  layer  for 
each  problem  encountered  in  this  study  have  been  transformed  by? mathematical 
techniques  into  ordinary  differential  equations.  Among  the  vast  number  of 
numerical  methods  available  for  solving  boundary-value  problems  of  ordinary 
differential  equations  the  Runge-Kutta  method  is  probably  used  most  frequently. 
Standard  Runge-Kutta  computer  programs  are  available  and  are  quite  easy  to 
utilize. 


. The  ordinary  differential  equations  of  Section  4.8  (Eqs.  4.66  and  4.67) 
for  the  equilibrium  flow  problem  can  be  expressed  in  terms  of  five  first-order 
ordinary  differential  equations  with  three  initial  conditions:  fs(0)  * 0, 

fs(0)  * 0,  g(0)  * gy,  and  two  boundary  conditions:  fgC0®)  * 1,  g(»y  = 1.  For 
the  frozen  and  nonequilibrium  flow  cases,  the  basic  equations  (Eqs.  4.69  to 
4.71)  can  be  written  as  seven  first-order  ordinary  differential  equations 
with  four  initial,  conditions:  fB(0)  « 0,  fB(0)  « 0,  g(0)  » gw  and  either 

z(0)  » constant  or  z'(0)  * constant,  and  three  boundary  conditions:  fg(w)  * 1, 
g(co)  = i and  z(oo)  = 1.  A procedure  is  now  given  for  finding  the  other  initial 
conditions,  fg(0) , g'(0)  and  z'(0)  (or  z(0)),  for  the  frozen  and  nonequilibrium 
cases.  A similar  procedure  can  be  used  for  the  equilibrium- flow  case. 


For  'he  calculations  we  assume  the  initial  conditions, 

f£(0)  - a 

g‘(0)  = b (Al) 

z*(0)  or  E(0)  = c 

and  define  the  following  equations, 

Fi(a,b,c)  = 1 - fg(») 

FE(a,b,c)  = 1 - g(»)  (A2) 

F3(a,b,c)  * 1 - z(») 

The  solutions  to  these  equations 

Fi(a,b,c)  - 0 

Fg(a,b,c)  » 0 (A3) 

Fa(a,b,c)  * 0 

can  often  be  found  by  a simple  procedure  known  as  the  Newton-Raphson  method. 

The  recurrence  relation  can  be  derived  from  the  Taylor  series  expansion  for 

£(xx) 

li%)  - f(x0)  + (xx  - xo)  f'(xx)  + ...  (A3) 

where,  f represents  Fj.,  Fg,  F3  and  x for  a,b,c. 


Al 


i 


From  Eq.  A3,  an  iterative  sequence  can  now  be  obtained  as  follows: 


ai+l  » a.  - VA 
Vi  * bi  ■ 

ci+l  " ci 


i 

£ 

x 

dFx 

ss: 

A 

= 

dFg 

c3a 

i 

8fe 

a? 

& 

i 

8f3 

*1 

Sf3 

cib^ 

dF3 

Fi(a.  ,b.  ,c. ) 

' 1*  1*  3/ 

SFx 

*T 

3a 

K 

ii 

Fe(a 

i’bi,0F 

Sf£ 

** 

dFe 

*1 

F3(a 

l’Vci> 

C>F3 

8f3 

<5c'i 

|a 

% 

Fi(ai, 

bi,ci> 

5fi 
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5f2 

<3aT 

X 

Fg^, 

bi’cF 

5f£ 

3f3 

i 

Fs^, 

bi>ci) 

Sf3 

i 

SFi 

mnB 

oa4 

i 

S>Fi 

i 

>bi'oi> 

n 

o 

O 

dFg 

dFg 

K 

Fs(ai 

>^i >ci) 

Sf3 

b?3 

Fs^ 

’VCi) 

A2 


m 


(A5) 


(AS) 


(A7) 


(A8) 


and  the  partial  derivatives  SFj/cJ&i  can  be  approximated  by 


+ h>Vci)  ‘ Fl^ai 


(A9) 


where,  h is  some  appropriate  step  size.  Similarly,  the  other  partial  derivatives 
can  be  evaluated. 


The  above  procedure  is  iterative,  since  an  initial  approximation  a = a0, 
b = bQ,  c = c0  to  the  root  of  Eq.  A3  has  to  be  made,  and  a sequence  of  approxima- 
tions bj_+^  and  Cj_+^_  are  generated  as  the  true  root  is  approached. 

The  initial  conditions  (Eq.  Al)  used  in  the  calculations  of  the  exact 
similar  solutions  of  Section  4.8  have  been  checked,  and  they  agree  with  the 
values  given  in  Refs.  23,  2k  and  63. 


■5 

j3 


% 

i 


A3 


,^-ji  -f 


